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FOREWORD 


The purpose of this report is to describe a recently completed 
program to design and manufacture an experimental transonic fan 
model featuring novel methods for noise reduction at the source. 

The program was conducted between 197^ - 1976 under contract 
NAS 3 -I 8512 issued by NASA Lewis Research Center, with Bolt 
Beranek and Newman Inc. (BBN) as the prime contractor and AVCO 
Lycoming as a subcontractor. The contract resulted from a NASA 
request for proposals (RPP) concerning CTOL aircraft engine fan 
source reduction concepts. The intent of the RFP was to identify 
advanced design concepts for reducing both rotor and stator 
sources which could be implemented with existing aerodynamic 
and structural design capabilities. The RPP encouraged proposals 
to reduce noise from high speed single stage fans. 

BBN proposed the use of "subsonic leading edge" rotor blades 
and variably swept stator vanes as the concepts to be investigated. 
The study and engineering work culminated in the fabrication 
of a 20-inch diameter fan stage to be tested for acoustic and 
aerodynamic performance at the Lewis Research Center, National 
Aeronautics and Space Administration, Cleveland, Ohio. 

Bolt Beranek and Newman Inc. (BBN), Cambridge, Massachusetts, 
served as the prime contractor with overall program responsibllty , 
as well as prime technical responsibility for the fan acoustic 
design, and other areas. The Lycoming Division of AVCO Corporation, 
Stratford, Connecticut, was a major subcontractor to BBN, with 
responsibilities in aerodynamic and mechanical design, and manu- 
facture of the fan hardware. Rotor blades and stator vanes were 
manufactured under subcontract by New England Aircraft Products, 
Farmington, Connecticut. 

Also included in the program were efforts to develop a 3- 
dlmenslonal compressible flow computer program to analyze the 
flow through the rotor, especially in the vicinity of the leading 
edge, and the investigation of the feasibility of using porous 
trailing edges on the stators to reduce broadband noise. The 
3-D flow program was discontinued at the time the rotor design 
was finalized, and the porous edge concept was not used because 
of the difficulties perceived in manufacture of small vanes from 
^^vallable porous metal materials. 

Numerous individuals at BBN and AVCO made significant contri- 
butions to this project. Mr. Richard Hayden served as project 


ill 



manager, and contributed to the acoustic design of the fan as well 
as other areas. Dr. Donald Bliss served as an associate project 
manager and had responsibility for the concept of the rotor blade, 
the rotor acoustic design, and the coordination of the aerodynamic 
design with AVCO. Mr. Bruce Murray also served as an associate 
project manager and supervised the mechanical design and manu- 
facturing aspects of the fan. The stator acoustic design was 
carried out by Dr. K.L. Chandiramanl, and Mr. Joseph Smullin. 

Drs. John McElman and John 0’ Callahan performed finite element 
stress analysis of the rotor blades, and Dr. O’Callahan contri- 
buted to numerical fluid mechanical analysis of the rotor flow 
field. 

At AVCO Lycoming, Mr. Pierre Schwaar served as the principal 
investigator and has primary program responsibility for the fan 
aerothermodynamical design, and for implementing the subsonic 
rotor leading edge concept and the acoustic design of the stator 
blades within operational structural constraints. Mr. Herbert 
Kaehler led AVCO’s work on structural analysis and Mr. John Banks 
supervised mechanical design and manufacturing activities there. 

Mr. James G. Lucas of the NASA Lewis Research Center’s V/STOL 
and Noise Division was the NASA Program Manager, and contributed 
valuable assistance in the mechanical design and manufacturing 
areas, and in the integration of the fan into NASA’s test facilitie 

This report has been designated as Bolt Beranek and Newman 
(BBN) Report No. 3332. 
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SUMMARY 


On current generation high bypass ratio turbofan engines, 
the fan is a predominant noise source which must be controlled 
to meet future aircraft noise goals. Of the various approaches 
to turbofan engine noise reduction, the most attractive is re- 
ducing the strength of the noise-producing elements at the 
source, thus avoiding weight and performance penalties associat- 
ed with various sound suppression approaches. 

In modern high bypass ratio turbofans, the fan thrust is 
achieved in a single fan stage, which usually requires super- 
sonic tip speeds of the fan rotor to produce the necessary 
pressure rise. In such fans, the predominant sources of noise 
are shocks radiated from the supersonically-moving rotor blades 
(called multiple-pure-tone [MPT] noise), and tones radiated from 
the rotor wake Interaction with stator vanes. 

In this program, two advanced noise reduction concepts were 
applied to the design of a 1.6 pressure ratio single stage fan. 
The goal of the design was to reduce the following acoustic 
sources: multiple pure tone noise, rotor-wake/stator-blade 

Interaction noise, and noise due to operating the rotor in dis- 
torted or turbulent Inflow. Unique nonradlal blading of the 
rotor and stator was used to achieve these goals. The rotor 
blade leading edges were swept so that the normal component of 
flow to the edge is subsonio at all points along the blade span, 
thus preventing the occurrence of leading edge shockwaves. The 
stator vanes were designed to minimize noise generated by rotor 
wakes incident on the blades by progressivly sweeping the vanes 
from root to tip in order to produce subsonic trace speeds for 
the unsteady loads along the span. Special aerodynamic and 
structural design considerations were required to assure the 
performance and integrity of this unusual blade and vane design. 

The rotor design using a blade concept with shock-free 
leading edges (except at points of inflection where weak conical 
shocks occur) is highly flexible in that a large family of 
blade shapes and leading edge contours may, in general, be 
used to achieve the noise reduction goal. The swept rotor de- 
sign is also attractive since it should perform equally well 
at off-design conditions if it has been designed to perform 
properly at the highest envisioned rotor speed. The swept 
edges also are compatible with reducing noise generation due 
to inflow distortion. 
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In the final design of the particular rotor ultimately 
constructed, a reversal of the sweep direction was required 
near mid-span to minimize stress levels in the blade. Once 
this was done, aeroacoustical-structural design iterations led 
to a blade with acceptable stress levels and no additional 
compromise in acoustic performance beyond the expected weak 
conical shock at the sweep reversal point. 

The aerodynamics of subsonic leading edge rotor cascades 
with supersonic absolute inflow velocities are not well known, 
and will clearly require further study. 

The concept of forcing the trace speeds of moving load 
distributions on stator vanes to be subsonic was introduced for 
the first time in this program. The design of a stator which 
uses this concept requires a controlled rate of axial sweep- 
back (or circumferential skew), the details of which depend 
heavily on the rotor wake field which varies with distance 
from the rotor. The selection of a stator vane number for a 
given rotor design is done with the familiar cutoff condition 
in mind; however, supersonic rotor tip speeds make it impossible 
to completely cut off the radiation at the tips of the stator 
vanes. No serious aerodynamic or structural problems were asso- 
ciated with the swept stator. The stator acoustic design pro- 
cedure is now well-defined in terms of flow parameters needed 
as inputs, but the ability to predict the necessary flow para- 
meters of the rotor wake field is presently limited. 



SECTION 1 


INTRODUCTION 


With the advent of high bypass ratio turbofan engines, 
and the associated decrease in exhaust velocity, the fan stage 
has become the dominant aircraft engine noise source. Therefore, 
fan noise reduction is a problem of primary importance in the 
ongoing effort to evolve quieter aircraft. Furthermore, it is 
increasingly important that any penalty in operating efficiency 
incurred by noise reduction methods be minimized. 

In general, noise reduction can be achieved in two ways: 

(1) reduction through the attenuation of propagating sound 
fields; and (2) reduction of the strength of the noise sources 
themselves. The first approach typically involves the use of 
absorptive duct liners and splitters, and possibly basic 
modifications to the inlet duct geometry. Because add-on 
features are required, and the duct length may be increased, 
the penalties associated with this approach are added weight 
and some direct reduction in aerodynamic efficiency. Further- 
more, there may be a degree of noise generation associated 
with some treatment modifications, such as in-duct splitters, 
particularly if the Inflow to the fan is disturbed. 

The second approach, which is the reduction of noise at 
the source, can be pursued in many ways. The basic fan design 
parameters can be chosen to give more favorable acoustic 
behavior. For Instance, the tip speed can be reduced, the 
spacing between the rotor and stator can be Increased, and the 
number of blades and vanes can be altered to prevent the 
propagation of certain duct modes. Whether these options 
can be exercised in a given case depends on the design con- 
straints on the performance and size of the system. 

Because, in most circumstances, acoustic considerations 
cannot dictate the choice of basic fan design parameters, other 
means of noise source reduction are worthy of consideration. 

These other means of source reduction necessarily Involve 
changes in the aerodynamic design of the blades and vanes. 

The design changes may occur either within the framework of 
conventional design practice, such as the use of optimized 
blade section properties, or may Involve the exploration of 
novel concepts. Although development of all the design data 
needed for Implementing novel concepts for noise source re- 
duction may be initially difficult, the noise reduction po- 
tential of a successful concept may greatly exceed the re- 
duction obtained by more conventional means. Of course, the 
final test of an acoustically successful concept must always 
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be whether any associated penalties in performance, complexity, 
and system integration can be overcome or, at least, justified 
in relation to the benefits. 

The subsonic leading edge rotor is implemented by tailoring 
(sweeping) the rotor leading edge to the mean inflow such that 
subsonic Mach number flow is achieved normal to the leading edge 
along the entire span, thus preventing shock generation. Previous 
use of partially-swept transonic rotors was done in an effort 
to reduce transonic drag rise and thus improve stage efficiency. 
Swept stators have been previously used to reduce noise, but the 
design concept implemented here involves tailoring the leading 
edge shape to a detailed estimate of the rotor wake field inci- 
dent upon the stator. 

The remainder of this report is organized to describe in 
detail the rationale for selection of the particular concepts 
(Sections 2 and 3)j details of the design procedure used on 
the swept rotor blades (Section 5) and stator vanes (Section 6), 
residual noise sources (Section 7), and facility integration 
(Section 8). Appendices contain a listing of aerothermodynamic 
design parameters (App. A), a discussion of geometric considerations 
for subsonic leading edge rotor blades (App. B), a detailed 
discussion of acoustical considerations in the stator design (App. 
C), discussion of empirical estimates of fan noise levels (App. D), 
and a useful algorithm for estimating trace speeds of rotor wakes 
on stator vane leading edges (App. E). 
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SECTION 2 

TRANSONIC FAN NOISE SOURCES 


This section summarizes the major noise sources and 
mechanisms encounted with transonic fans. Typical design 
characteristics of single stage transonic fans are summarized 
in Table 1. 

TABLE 1. TYPICAL CHARACTERISTICS OF SINGLE STAGE TRANSONIC 
FANS. 

Pressure Ratio Range 1.4 - 1.8 

Tip Speed 300 to 600 m/s (1000-2000 ft/sec) 

Relative Rotor Tip Mach No. 1.1 - 1.8 
Rotor Inlet Hub/Tip Ratio .35 - .50 
Stator Hub Mach No. .8 

The most important noise sources, which Involve both the 
rotor and stator, are : 

Shockwave noise from the supersonic portion of the rotor 
blades, often called multiple pure tone (MPT) noise. 

Rotor/stator Interaction noise caused by unsteady loading 
due to aerodynamic Interaction (tonal and broadband noise). 

Noise caused by unsteady loading on rotor blades inter- 
acting with Inflow distortions and turbulence (tonal and 
broadband noise). 

A brief elaboration on each of these sources is now provided. 

2.1 Shockwave Noise 

When the relative flow past the rotor becomes supersonic, 
the propagation of shock waves out of the inlet duct becomes an 
important noise source. The upstream propagation of waves from 
blade rows with detached and attached shock wave patterns is 
shown in Fig. l,,(from Ref. 1). Because the pressure field must 
satisfy a periodicity condition, expansion waves occur in the 
regions between the shock waves. 

Several Investigations (Refs. 2 through 7) have shown that 
nonlinear effects are an Important factor in the upstream shock 
propagation process. Because nonlinear attenuation occurs more 
rapidly for higher initial levels, an increase or reduction of the 
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Expansion 

Wave Region Oblique Shock 



(a) Detached Shock Wave Pattern 



Rotor Motion Swallowed Shock 


(b) Attached Shock Wave Pattern 

FIG. 1. POSSIBLE SHOCK WAVE CONFIGURATIONS FOR 
ROTORS IN SUPERSONIC FLOW. 
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shock strength at the blades does not produce a comparable 
increase or reduction of levels at the end of the inlet duct, 
or in the far field. This effect is strongest when the wave 
train in the duct is well ordered and can be considered nearly 
one-dimensional in character. The Important consequence of 
this effect is that very substantial levels of source reduction 
must be achieved to guarantee a worthwhile reduction in level 
in the far field. 

Another Important consequence of nonlinear propagation is 
the redistribution of the shock noise spectrum from blade 
passage frequency and its harmonics to the rotor shaft rota- 
tion frequency and its harmonics. This redistribution occurs 
because of blade-to-blade differences in the initial strength 
and position of the shock waves. These blade-to-blade differ- 
ences are caused by variations in manufacturing tolerances 
that may affect the circumferential location, setting angle, 
thickness, and camber of the blades. Because the shock train 
structure is Inherently unstable to perturbations in strength 
and position, these Initial disturbances need not be large. 

As an example, when periodic variations in shock strength 
occur, the stronger shocks tend to overtake and dominate the 
weaker shocks because of nonlinear effects. Because the 
variations in strength are caused by blade-to-blade differences, 
they are periodic in the shaft rotation speed. Thus, as the 
wave train propagates, the harmonics of shaft speed become 
increasingly Important relative to the harmonics of blade 
passage frequency. Pig. 2 shows the redistribution of energy 
from blade passage frequency to shaft rotation frequency as 
the result of an initial amplitude perturbation to one shock 
in a wave train. Figure 3 shows sketches of typical noise 
spectra for a subsonic fan, which has no shock noise, and for a 
supersonic fan, where the tones at the harmonics of shaft speed 
are clearly present. Clearly the multiple pure tone noise due 
to shock wave propagation Is a major noise problem. 

2.2 Rotor/Stator Interaction Noise 

Unsteady aerodynamic loads on rotor blades or stator vanes 
produced by the aerodynamic interaction between the rotor and 
stator are an Important source of both tonal and broadband noise. 
The main causes of the aerodynamic interaction are the Interference 
with the potential flow pressure and velocity fields and the 
interaction with the viscous and turbulent wakes from upstream 
blades. The potential field Interaction that produces tonal noise 
at the harmonics of the blade passage frequency can be virtually 
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SOUND PRESSURE LEVEL SOUND PRESSURE LEVEL 



a ) TYPICAL NOISE SPECTRUM FOR A SUBSONIC 
TIP SPEED FAN 



HARMONICS) 


FREQUENCY 

b) TYPICAL NOISE SPECTRUM FOR A SUPERSONIC 
TIP SPEED FAN 


FI6. 3. TYPICAL FAN NOISE SPECTRA FOR SUBSONIC 
AND SUPERSONIC TIP SPEEDS. 
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eliminated by providing adequate spacing between the rotor and 
stator. Increasing the spacing on a high by-pass ratio fan 
stage is usually practical and does not involve a severe aero- 
dynamic penalty. The interaction of the stator vanes with the 
•'mean component" (steady velocity deficit) of the rotor wakes 
produces tonal noise at the harmonics of blade passage fre- 
quency, while the interaction with the wake turbulence produces 
broadband noise. Increasing the spacing between the rotor and 
stator also reduces - but does not necessarily eliminate - this 
noise source. 

2.3 Inflow Distortion Noise 

The inflow to the fan rotor typically exhibits a degree 
of spatial nonuniformity and a certain amount of turbulence. 
Sound is produced by unsteady loads on the rotor blades 
operating in this disturbed inflow. Steady spatial nonuni- 
formity causes tonal noise to be produced at the harmonics of 
blade passage frequency, and the presence of turbulence 
produces broadband noise. However, if the turbulence scales 
are sufficiently long in the streamwlse direction, then many 
blades will interact with a given disturbance in a similar 
manner, producing peaks in the noise spectrum at the harmonics 
of blade passage frequency. Because the basis for this noise 
source is a random process, the amplitude of these peaks will 
vary in time in a random manner. Inflow distortions have been 
shown to be a potentially important noise source in static fan 
test facilities. Their importance in an actual flight enviro- 
nment is less certain, since the effect of forward motion is 
usually to reduce certain types of inflow distortion. 
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SECTION ^ 

NOISE SOURCE REDUCTION CONCEPTS 


In this section, the concepts for the reduction of rotor 
and stator noise sources are described. A review of the de- 
tailed analysis and design procedures associated with the Im- 
plementation of these concepts In the present program Is post- 
poned to the sections later In the report dealing with detailed 
design. 

3.1 Rotor Noise Reduction 

As discussed In the previous section, two noise sources 
associated with the rotor are multiple pure tone noise due to 
shock waves and Inflow distortion noise. This section describes 
a concept which has the potential to substantially reduce multi- 
ple pure tone noise. As an additional advantage, this concept 
will also help reduce the problem of inflow distortion noise. 

In principle, upstream-propagating shockwave noise can be 
reduced by designing for careful alignment of the relative 
velocity, w, with the suction surface near the rotor blade lead- 
ing edge, as shown In Pig. 4a. However, completely shockfree 
entry Into the blade row cannot be achieved in conventional 
blading because of the finite thickness of the blade leading 
edge. The effect of thickness is Illustrated in Fig. 4b. More- 
over, since the relative Inflow direction varies with the opera- 
ting conditions, the proper alignment cannot be maintained In 
off-design operation, nor in the presence of Inflow distortions. 
Thus, this concept presents several practical difficulties for 
application to aircraft fans which do not operate at a single 
design point. 

A different approach to obtain shockfree entry Into a 
blade row Is now described. It Is believed that this approach 
does not suffer from the shortcoming of the more conventional 
approach just described. Consider a blade whose leading edge 
is swept relative to the local Inflow velcocity vector. The 
leading edge would in general appear swept when viewed from 
the side and skewed when viewed from the front. If the leading 
edge is swept such that the Mach number of the relative flow 
component normal to the leading edge is everywhere subsonic, 
a shockless leading edge results. In wing theory, this Is 
referred to as a "subsonic leading edge in supersonic flow" 

(See, for Instance, Ref. 8). In rotating applications, the 
radial variation In relative Mach number makes it possible, 
in principle, to completely avoid upstream shock wave propagation 
by using leading edge and surface generating line sweep which 
varies from hub to tip. In practice, structural constrains force 
some design compromises. In the present design, the structural 


9 



Upstream 


FIG. 



I 


c 
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compromise entails the presence of a train of conical shocks 
upstream of the rotor associated with a sweep discontinuity 
in the leading edge. From the standpoint of preventing shock 
noise, the design can be made insensitive to operating condi- 
tions, relative flow alignment, and inlet distortions by de- 
signing the sweep distribution for the highest relative Inflow 
Mach number to be expected; thus ensuring a lower subsonic 
normal Mach number component for off-design conditions. This 
Insensitivity is considered to be a major asset. 

The underlying aerodynamic idea is now reviewed. Figure 5a 
shows a swept wing of infinite extent subject to an Incident 
supersonic flow. Since there is no spanwlse variation in the 
wing geometry, the axial component has no effect. The aero- 
dynamic forces are determined entirely by the component of the 
flow normal to the wingspan. If the component normal to the span 
is subsonic, then there are no shock waves associated with the 
flow over this wing. Of course, to be completely shockless, the 
normal component must be sufficiently subsonic that transonic 
flow effects do not occur in the normal flow plane. The only 
effect of the axial component is in the structure of the viscous 
boundary layer on the wing surface, but this is not related to the 
presence or absence of shock waves. The same ideas are appli- 
cable, of course, to an infinite span sweptback cascade. 

Fig. 5b shows a finite span wing sweptback to have subsonic 
leading edges. The aerodynamic behavior is now considerably 
more complicated. In particular, the presence of conical 
shocks at the front and rear of the wing root and at the 
rear of the tips is unavoidable. These Isolated points on the 
wing are discontinuities in the otherwise subsonic edges. 

The conical shocks are, however, weaker than their two- 
dimensional counterparts and, because of their three-dimensional 
nature, decrease in strength with distance from their point 
of origin. 

The application of a subsonic leading edge to a fan blade 
is Illustrated in Fig. 5c. This illustration is simplified to 
its essential form, showing only the radial change in Mach 
number. The actual process is nonplanar because of the change 
in direction of the Inflow with radial location. The particular 
case illustrated applies to a transonic fan, since part of the 
Incident flow is subsonic . Then the leading edge can be made 
completely shockless even though the blade is of finite extent 
assuming that one is able to predict and accommodate the effects 
of spanwlse flows (Ref. 9 ). 
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The local leading edge sweep at each radial station is chosen 
to be greater than the Mach angle of the local flow, l.e., 
the swept edge must lie within the local Mach cone. This assumes 
that the normal flow to the leading edge is everywhere subsonic. 
Because of the gradient in Mach number, the Incident flow is sub- 
sonic at the base of the blade so a shock cannot emanate from 
this point (unlike the wing root in Pig. 3b). Hence, the blade 
leading edge can be entirely shockless, except for the effects of 
aerodynamic Interference between the blade tip and the shroud 
which produce conical shocks . If the fan were completely super- 
sonic, a conical shock should also occur at the root of the blade. 
By designing the leading edge and the other generating lines of 
the forward portion of the blade surface to be subsonic for the 
situation that produces maximum relative flow Mach number, the 
edge will remain subsonic under all other operating conditions. 

The blade leading edge would usually be designed to have a con- 
stant normal velocity (Mach number) component at all points along 
the span at radii (from the hub) greater than that at which the 
critical normal Mach number, is reached. The critical Mach 

number is that normal Mach number (<1) at which thickness ef- 
fects would cause the flow to become transonic. 


In addition to sweepback. Pigs. 6a, 6b, and 6c show swept 
forward and compound sweep blades that are also possible 
configurations. All of the blade configurations must have a 
conical shock at the tips caused by aerodynamic interference 
with the shroud. The compound sweep blade will also have a weak 
conical shock at the discontinuity in sweep, which is positioned 
somewhere along the leading edge (assuming the discontinuity 
lies in the region of supersonic relative Inflow). Although 
the compound sweep blade has the acoustic penalty of introducing 
a weak conical shock, it offers other definite advantages. 
Structural considerations provide the most severe constraint 
to the design of high speed fans with swept blades. Palrly 
large excursions of the leading edge are required to implement 
this concept. It should be noted that the family of three- 
dimensional curves that satisfies the subsonic leading edge 
condition is not unique and therefore considerable latitude exists 
to determine structurally optimum shapes. Pigure 7 shows the 
type of conical shock wave pattern for a compound sweep blade. 

The blade in the sketch closely resembles the design developed 
during the course of the project being described. 

Pigure 8 compares the operation of a moderately loaded 
blade row with and without subsonic leading edges in supersonic 
flow. As explained above, the subsonic edge region allows 
shock-free entry into the blade row. The blade rows are 
Identical except for the addition of a subsonic leading edge 
region in one case. The front surface of the blade must be 
designed so that any shocks generated on the suction surface of 
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FIG. 6 



B. Swept Forward 



C. Compound Sweep 


. FRONT VIEW OF SOME POSSIBLE BLADE CONFIGURATIONS 
WITH SUBSONIC LEADING EDGES. 
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FIG. 7. CONICAL SHOCK FIELD FROM A ROTOR BLADE 
WITH A COMPOUND SWEEP LEADING EDGE. 


15 



(a) OPERATION OF A CONVENTIONAL MODERATELY LOADED 
BLADE ROW 



(b) OPERATION OF A BLADE ROW WITH SWEPT (“SUBSONIC") 


LEADING EDGES 


FIG. 8. COMPARISON OF THE OPERATION OF A MODERATELY LOADED 
BLADE ROW WITH AND WITHOUT SUBSONIC LEADING EDGES. 
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the blade are formed sufficiently far back that the disturbance 
is entirely contained in the blade row, even during off-design 
operation. 

Using a swept leading edge also helps reduce the response 
of the rotor to inflow distortions, because the magnitude of 
the response is largely determined by the velocity component 
normal to the leading edge (Ref. 10). The effect of inflow 
distortion is most Important near the tip of the rotor where 
the relative velocity is highest. Fortunately, the concept 
for sweeping the blades requires the most sweep near the tip. 

3.2 Stator Noise Reduction by Leading Edge Sweeping and Blade/ 

Vane Number Selection 

Although increasing the spacing between the rotor and stator 
leads to some noise reduction, the aerodynamic interaction be- 
tween the rotor wakes and stator vanes remains an Important noise 
source. Further reduction by conventional means can be achieved 
by choosing the proper number of blades and vanes to cut off 
many of the acoustic spinning modes in the duct (Ref ..ii'» . When 
the rotor tip speed is subsonic, the blade and vane numbers can 
be chosen so that all the spinning modes at blade passage fre- 
quency, and at least some of the modes at higher harmonics, are 
cut off. However, if the rotor tip speed is supersonic, at least 
one spinning mode at blade passage frequency cannot be cut off, 
regardless of the choice of blade and vane numbers. Since super- 
sonic spinning speeds often occur on transonic fan designs, other 
means of stator noise reduction are of considerable Interest. 

Figure 9a Illustrates the Interaction of a row of stator 
vanes with rotor wakes when viewed on a surface of constant 
radius from the fan axis. The wakes can be described as flow 
regions with an average velocity W lower than the velocity of 
the adjacent fluid, upon which a turbulent perturbation velocity 
field Aw is superimposed. 

Figure 9b shows a sketch of a three-dimensional wake/vane 
Interaction in a fan. The structure of the viscous, usually 
turbulent, wakes that trail each rotor blade is complex. How- 
ever, on the average, these wakes can be considered as being 
convected with the mean flow in which they are Imbedded. The 
nature of the downstream mean flow is such that the convection 
process will distort the wakes from their original shape; namely, 
the downstream flow is distorted both axially and circumferentially 
across a given radial path, leaving the downstream pattern of the 
wake disturbance very much altered from the pattern at the rotor 
trailing edge. Suppose, for instance, the rotor is designed to 
give a mean flow that has a uniform axial velocity distribution 
and a free vortex tangential velocity distribution. Assuming 
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(a) The Interaction of the Stator Vane Row with the Mean and Unsteady Rotor 
Wake Components as Seen on a Constant Radius Surface. 



(b) A Sketch Showing the Three-Dimensional Nature of the Rotor-Wake/ 
Stator- Vane Interaction. 
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FIG. 9. THE CHARACTERISTICS OF ROTOR-WAKE/STATOR-VANE 
INTERACTION. 








the wakes are radial at the rotor trailing edge, it is clear 
that the tangential velocity component will act to skew the 
wakes over, with the hub region ' leading the tip region. This 
situation is Illustrated in Pig. 9b. In this case, the inter- 
action of a given wake with a given stator vane does not occur 
simultaneously all along the stator vane span. Instead, the 
Instantaneous spanwlse interaction region of a single rotor 
wake will extend over only a portion of any one vane and will 
sweep along the vane leading edge, beginning at the hub and 
ending at the tip. The skewing of a wake due to convection 
by the downstream mean flow can be sufficient to involve 
simultaneously portions of several stator vanes. 

The shape of wake and the magnitude of its velocity components 
vary from hub to tip. To complete this picture of the downstream 
flow field, one must consider the unsteady velocity components 
which account for the turbulent structure of the wakes and for 
any other sources of Inhomogenietles in the flow, e.g., inlet 
flow distortions, large-scale flow instabilities, and blading 
errors. In general, the statistical properties of these unsteady 
components can be expected to vary axially, circumferentially, 
and radially. 

Both the mean and unsteady velocity components of the wake 
flow induce unsteady loads on the stator vanes. The mean com- 
ponent will produce a load distribution that travels from hub 
to tip, changing shape and amplitude in accordance with the 
radial variation of the mean flow properties and wake strength, 
width, and skew. Imposed on this traveling load distribution 
will be the unsteady effect of the turbulent structure of the 
wake. The end result of all sources of unsteady loading on 
the stator vanes is to produce tonal and broadband noise. The 
tonal noise is usually considered to be the more Important 
noise source. The speed at which the point of interaction of the 
flow disturbance with the vane travels along the span is called 
the trace speed. 

A particular source of unsteady loading will produce no 
significant acoustic radiation if it satisfies a subsonic and 
non-accelerating trace speed criterion along the vane span. The 
trace speed concept has been previously recognized for the problem 
of helicopter-blade/vortex interaction by Widnall (Ref. 12) although 
it has not been generally recognized in the study of fan noise. 

The interaction of the wake with the vane produces a load 
distribution that travels along the vane. Suppose the vane is 
much longer than an acoustic wavelength. Following the trace 
of a phase front of this load distribution, acoustic radiation 
can occur along the vane span if the magnitude of the load changes. 
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the phase speed changes with time, or if the phase speed Is 
supersonic. For Instance, in fan noise analyses the rotor-wake/ 
stator-vane interaction is usually assumed to be two-dimensional 
(corresponding to infinite spanwise trace speed). The conditions 
mentioned above are necessary for radiation but not sufficient. 

The interaction with the acoustic field produced by the other vanes 
must also be considered before the actual occurrence of acoustic 
radiation can be established. Therefore, regions along the 
stator vane span can be expected to be poor radiators if the 
phase speeds are subsonic, nearly constant, and local levels do 
not vary rapidly. Other regions may or may not be efficient 
radiators depending on the behavior of the distribution of 
sources elsewhere on the stator. Furthermore, end effects at 
the hub and tip (within approximately one half an acoustic wave- 
length of the ends) makes these regions potential radiators. 

These considerations are discussed in Appendix C, and justified 
in detail in Bliss, et al . , (Ref. 13). 

Understanding the rotor-wake/stator- vane Interaction and 
the criteria for radiation from the span of a single vane 
suggests ways in which the vane configuration can be altered 
to achieve noise reduction. The vane should be shaped so that 
loads traveling along the span move at a constant subsonic 
speed. Assuming that the amplitude of the load distribution, 
moving with a phase front, is essentially constant, then 
radiation from the vane span will not occur (except for end- 
effects). The condition of a constant subsonic spanwise trace 
speed can be achieved by sweeping or skewing the stator vanes, 
as Illustrated in Fig. 10. In this illustration, the lines of 
constant phase can be considered to be the Intersection of the 
rotor wakes with the olane of observation (e.g., the r-9 plane 
in Fig. 10a, and the r-z plane in Fig. 10b). Except for 
the effect of shape changes, these lines travel at constant 
speed (rotational in the r-6 plane and rectilinear in the r-z 
plane) because of the rotation of the rotor. The speed at 
which a phase front traces the leading edge of the stator vane 
depends on the shape of the leading edge and the shape of the 
phase front. Clearly the trace speed can be controlled by 
either sweeping or skewing the stator vane. With this approach, 
radiation from the stator span can be prevented, leaving only 
acoustic radiation from end effects at the hub and tip of the 
vane. Radiation from the hub region can be cut off by the 
proper choice of blade and vane numbers, provided that the 
rotation speed of wakes at the hub is subsonic. Since the 
rotation speed of wakes at the stator tip will usually be 
supersonic for a transonic fan, the radiation from tip end 
effects can never be entirely cut off. Note that the rotor 
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Skewed Stator Vane 


FIG. 



(a) Schematic Axial View (r-0 plane)-Skewed Vane 


Lines of Constant 
Phase in Rotor Wake 



(b) Longitudinal Section (r-z Plane)-Swept Vane 

10. SCHEMATIC OF PHASE SHIFTS BETWEEN ROTOR WAKES 
AND SKEWED/SWEPT STATOR VANE. 
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wake pattern rotates with the same angular velocity as the 
rotor. Thus, at any given radius at any downstream location 
between the rotor and stator, the rotation speed of the wake 
pattern is simply, fir, which is different than the swirl velocity 
component. This can be best visualized from rotor fixed co- 
ordinates from which the wake pattern appears "frozen." 

Another, but related, way to view the effect of sweeping 
or skewing the stator vanes is as follows. Tyler and Sofrln 
(1962) have shown that for a given circumferential mode 
number, m, and hub-to-tip ratio, v, the radial structure of 
an acoustic spinning mode can be described by an infinite 
series of characteristic functions. The functions in this 
series differ according to their radial order, y, l.e., each 
function has a different number of nodes in the interval 
between the hub and tip. The spinning speed at which each of 
these functions begins to radiate is always supersonic and 
increases with increasing radial order. Therefore, at a given 
supersonic spinning speed and fixed m and a, only a certain 
number of the functions corresponding to the lowest radial order 
will not be cut off. Vanes can be skewed or swept so that 
the number of wakes on a given vane is increased, raising the 
radial order of the load distribution on the vanes. The 
acoustic energy is thereby redistributed to higher radial 
orders, some of which will be cut off. The relationship between 
duct mode cut off and the constant subsonic trace criterion 
is discussed by Bliss, et al. , (Ref. I3), and in Appendix C. 



SECTION 4 


FAN STAGE DESIGN SUMMARY 


An experimental transonic fan stage was designed and con- 
structed using the noise reduction concepts explained in the two 
preceding sections. The fan uses compound sweep rotor blades 
designed to have "subsonic leading edges" in the region of super- 
sonic relative Inflow. The stator vanes were swept back to 
achieve a constant subsonic trace speed of rotor wakes along the 
vane span. Figures 11a, b and c show photographs of the actual 
fan stage. A cross-sectional view of the fan as it will appear 
when installed in the test facility of NASA Lewis is shown in 
Pig. 12. As Indicated in the illustration, the fan will be 
tested in both forward and reverse Installation arrangements in 
order to measure both the fore and aft noise characteristics. 

The design data for the fan stage is summarized in Table 2. 

In the remainder of the report, the detailed design procedures 
used in the development of the fan stage are described. 

TABLE 2. FAN STAGE DESIGN SUMMARY 


Stage Characteristics: 

Stage Pressure Ratio, P^/Pi =1.6 

Mass Plow Rate, W = 31«2 kg/s ( 68.8 Ib/sec) 

Specific Mass Flow Rate : (referred to annular area at 

rotor inlet) 

W = 199 . 03 kg/s*m^ (40.76 Ib/sec* ft ^ ) 

0.0 

Polytropic Stage Efficiency, n =0.86 


Rotor : 


28 Compound Sweep Blades 

Leading Edge Normal Mach Number = 0.91 

Tip Speed = 480 m/s (1575 ft/sec) 

Relative Tip Inlet Mach Number = 1.593 
Rotor Inlet Tip Radius = 249 mm (9.o03 in) 
Rotor Inlet Hub-Tip Ratio = 0.442 
Rotor Pressure Ratio, P 2 /P 1 =1.64 

Stator : 

59 Swept Back Vanes 

Sweep Angle = 25^ At Root, 40 At Tip 
Stator Inlet Mach Number = 0.80 
Stator Pressure Loss AP^_^/Pj = .025 
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(a) 45° VIEW OF ROTOR AND STATOR ASSEMBLY FROM FRONT OF FAN 

FIG. 11. PHOTOGRAPHS OF THE EXPERIMENTAL FAN STAGE 










SECTION 5 

DETAILED ROTOR DESIGN 


This section deals with all aspects of the detailed 
analysis and design procedures associated with the fan rotor. 

The aerodynamic design of a transonic rotor having blades with 
"subsonic" leading edges differs significantly from conventional 
design practice because the acoustic, aerodynamic, and 
structural requirements interact strongly with each other from 
the very beginning of the preliminary design phase. Therefore 
it was necessary to conduct numerous aerodynamic-acoustic- 
structural design iterations to optimize and finalize a 
rotor configuration satisfying all the design requirements. 

The overall design point data for the rotor are listed in 
Table 2 of the previous section. 

5.1 Aerodynamic Design 

Certain differences are to be expected in the aerodynamic 
behavior of a rotor with subsonic leading edges. Since the 
entry into the blade row is shock free, any shocks that occur 
must remain within the blade row under all operating conditions 
because the edge region cannot support a shock system. Further- 
more, the effects of sweep may introduce other three-dimensional 
flow phenomena which are not present in a conventional blade 
design. Given these facts, the rotor aerodynamic design was 
undertaken using the best conventional design practice combined 
with an anticipation of the most important effects of swept 
edges. The design was carried out primarily with the use of 
an axisymmetrlc flow computer program. Conventional (two- 
dimensional) methods for analyzing the flow behavior within the 
blade row are not really adequate for the three-dimensional case 
of blades with swept "subsonic" edges. To handle this problem 
analytically requires a more general approach. Some work was 
done to adapt a new fully three-dimensional computer code to 
the analysis of flow through the blades with swept edges, but 
was discontinued due to schedule requirements. 

An important question in the design of a rotor with "sub- 
sonic" edges is related to its surge margin. Because the edge 
region cannot support a shock system it was felt that the 
surge margin might be reduced. Such a reduction would occur, 
if the effective rotor operating range were limited by the 
condition that the shock system remain within the covered 
cascade region. The flow configuration in which the shock 
system must remain within the covered cascade, however, does 
not yield the maximum static pressure rise achievable in a 
conventional transonic-supersonic rotor. Consequently, in a 
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stage where surge Is not triggered prematurely by the stator 
flow conditions, a rotor with subsonic leading edges might 
result in a decrease of the surge margin as compared to a 
conventional design. 

Since the rotor aerodynamic loading essentially depends 
upon the rotor static pressure rise, the selection of the 
meridional flow path was the main design step taken to achieve 
the desired loading levels. Meridional channel conicity and 
curvature through the rotor section were traded off in several 
preliminary design attempts. The flow calculations were per- 
formed by means of a code which solves the general equation 
of radial equilibrium on straight axial or slanted stations 
for the axlsymmetrlc flow case taking into account the radial 
variation of the blade efficiency. The polytropic efficiency 
n assumed for the rotor blading is shown in Pig. 13 j where n is 
derived from 



It was found that the comparatively large channel conicity across 
the rotor section required by the high design pressure ratio 
P^/Pj = 1.6 and the wall curvature needed at rotor exit to pre- 
vent excessive channel contraction in the free space between 
rotor and stator, combine to shift the maximum rotor static 
pressure rise from the tip towards the midspan location, where 
the shock system has the greatest tendency to move upstream 
into the uncovered cascade region because of the lower relative 
inlet Mach number. The main preliminary design effort con- 
sequently was directed toward minimizing the static pressure 
rise at the critical midspan location. 

The optimum channel configuration is shown in Fig. 14. 

The flow conditions are summarized on the Aero design program 
(R-121) input and output printout attached in Appendix A. 

Figure 15a shows the distribution of the rotor static 

pressure ratio ?^/? over the channel height, together with 

the relative inlet Mach number M and the corresponding nor- 

^ _ 

mal shock pressure ratio Pj/Pj, which is roughly equivalent to 
the static pressure ratio obtained in the front portion of a 
conventional cascade with a normal shock attached to the leading 
edge . 
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MERIDIONAL FLOW PATH. 




Three types of operating conditions can be distinguished 
along the rotor blade span. 

(a) From the hub to the sonic radius r ... the rotor 

i(M = 1)» 

1 

static pressure rise is achieved essentially by subsonic relative 
flow deceleration and centrifugation. 

(b) From r^^j^ _ to r^^p ^p ^ the rotor static 

2111 

pressure rise must be achieved through a channel-contained normal 
shock or a pseudoshock system followed by subsonic relative flow 
deceleration. The radial distribution of the rotor static pres- 
sure ratio P 2 /P 1 determined how far this operating condition ex- 
tends beyond the sonic radius. In the present case, it extends 
roughly from the \2% to the h0% mass flow streamline, or from 
20% to 53% of the span, l.e., slightly beyond the point of 
sweep reversal. The maximum inlet Mach number in this blade 
section remains below the I.3 level at which the interaction 
of normal shock with a turbulent boundary layer produces ex- 
tensive flow separation (P/P = 1.8). If minor flow separation 
does occur in the upper portion of this region, the flow will 
reattach to the blade because of the large solidity provided 
in the vicinity of the point of sweep reversal. Consequently, 
it is expected that the design flow conditions will be obtained 
over this critical span section by a shock configuration located 
in the forward, yet still covered portion, of the cascade. 

(c) In the upper blade section, P 2 /P 1 IS smaller than P^/Pj, 
and the shock system consequently moves progressively toward 
the rear portion of the cascade. Since no shock is attached 

to the leading edge, the flow conditions are essentially similar 
to those in the diverging section of a converging-diverging 
nozzle in the supersonic off-design operating range. 

Figure 15b schematically shows the meridional projection 
of the rotor blade and the anticipated shock/pseudoshock 
interception area on its pressure and suction sides. The main 
question pertains to the rotor surge margin, i.e. the extent 
to which the tip region will be allowed to increase its 
pressure ratio beyond the design value by forward shifting of 
the shock configuration before (i) flow separation occurs at 
the hub, or (li) the shock system at midspan is forced into 
the uncovered cascade region. 
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FIG. 15(a). ROTOR FLOW CONDITIONS: SPANWISE DISTRIBUTION 

OF STATIC PRESSURE RATIO AND INLET RELATIVE 
MACH NUMBER (Design Point). 



FIG. 15(b). ROTOR BLADE SHOCK/PSEUDOSHOCK INTERCEPTION 
AREA (Schematic). 
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5.2 Description of the Aero-Structural Design Interaction Problem 

The problem of achieving acceptable stress levels is much 
more difficult for a rotor with swept leading edges than for a 
conventional design. The aerodynamic and structural require- 
ments for the rotor blade are therefore closely coupled. Within 
the aerodynamic constraints, a number of design iterations were 
required to achieve acceptable stress levels and to optimize the 
design. The major aerodynamic constraints are that the rotor 
meet the design performance requirements and that the normal 
component of flow to the leading edge be maintained at a certain 
subsonic value. Because of the gradient of relative inflow Mach 
number, the angle of sweep must increase toward the tip to meet 
the condition of a subsonic normal component. In the present 
case, the maximum normal component Mach number was chosen to be 
0.91 along these leading edges (actually lower near the hub). 

The value of 0.9I was chosen as a goal since it represented a 
normal Mach number sufficiently below sonic to avoid thickness- 
related shocks. Lower values can be chosen, but the severity of 
the blade leading edge excursions increase as the normal Mach 
number is lowered. For the fan design tip speed, the excursions 
of the swept leading edge are large and it was necessary to use 
a compound sweep configuration to minimize bending stresses. The 
major variables available to control blade stresses are the loca- 
tion of the sweep reversal point, the local section properties of 
chord length, maximum thickness, thickness distribution, and the 
stacking of the blade sections. 

Because of the large leading edge excursions, the centers 
of gravity of the blade sections can no longer be stacked on a 
radial line. In addition to the centrifugal tensile stresses, 
large bending moments about both principal axes of inertia of 
the blade sections were found to occur (Fig. 16 ). Achievement 
of acceptable stress levels required the use of a carefully 
chosen sweep reversal point and the development of an effective 
nonradial stacking procedure. 

Typically, the most critical problem was the bending 
moments about the minor axis of inertia, and a special 
stacking procedure was used to minimize these moments. A 
near-optimum procedure for nonradial stacking is as follows. 

The blade sections were stacked starting at the tip and moving 
inward. The addition of each Incremental blade section was 
made so that the center of gravity of the entire portion of the 
blade above this section falls on the axis of minimum inertia 
of the new section. The center of gravity of the new upper 
portion was then reevaluated before the next incremental section 
was added in the same manner. This procedure nearly minimizes 
the critical bending stresses around the axis of minimum 
inertia. The result is not completely optimum because of the 
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complexity of the actual situation in which the stresses are 
determined by the complex interaction of many effects. Further 
improvements were made by iterative changes around the result 
of the above stacking procedure, particularly with the inten- 
tion of relieving local stress concentrations. To the extent 
that high stresses arise due to bending around the axis of 
maximum inertia, these can be relieved largely by changing the 
location of the sweep reversal point and varying the local 
section chord and thickness. 

5.3 Determination of the Subsonic Rotor Leading Edge Geometry 

At each leading edge point, the relative Mach vector 

1 

defines a Mach cone. To a prescribed value of the subsonic 

velocity component M perpendicular to the leading edge, 

w 

'L 

there corresponds a coaxial cone with smaller aperture. The 
subsonic leading edge elements must only satisfy the condition 
that they lie on such cones. Their direction otherwise is 
arbitrary . 

Referring to Fig. 17, a particular sweep direction can be 

defined by specifying that each leading edge element be swept 

in the plane formed by the relative inlet velocity, W, and the 
radius (W-r) plane. This yields the shortest leading edge line 
from hub to tip, since it maximizes the radial projection of 
every leading edge segment. 

Sweeping in the W-r plane however, does not result in a 
blade with minimum stresses. The resulting stacking of the 
CG’s of the profiles in fact was shown to generate substantial 
bending moments around their axis of minimum Inertia. The 
main parameter used to minimize bending stresses is the lateral 

sweep angle, V, between the radial plane passing through the 

leading edge element dl and the W-r plane. The situation is 
Illustrated in Fig. 17. The geometric analysis used for this 
design is described in Appendix B, and only some pertinent 
results are cited below. They are expressed by the two follow- 
ing equations for the cylindrical coordinates and Zj^ of the 

leading edge points in function of the relative flow angle 3^, 

the lateral sweep angle v, the slope e of the relative velocity 

w 

1 

and the projection y" of the Mach cone angle y in the W-r plane 
(see Appendix B) : 
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The relation between y" and y is given by the formula 
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In the above relations the (+) sign applies for backward, 
the (-) sign for forward sweep. 

The formulae define a sonic leading edge, l.e., leading 
edge points lying on the Mach cones of the adjacent points. 

A subsonic leading edge is simply obtained by using in the 
formulae y values corresponding to relative Mach numbers 


increased by a factor f = l/„ , l.e. M^ = M, 

W 


''■L 


w ,/M 


where 


w . 
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Is the subsonic Mach number of the relative velocity 

component perpendicular to the leading edge. This simple re- 
lationship is illustrated in Fig. l8. 

The second design parameter used to minimize the bending 
stresses was the sweep reversal radius. By proper selection 
of the point of sweep reversal, the center of gravity of the 
blade can be located in such as way as to project radially on, 
or near, the axis of maximum inertia of the hub section. From 
a structural viewpoint, the compound sweep blade of Fig. 6 
could be considered as a blade with hub and tip sections de- 
signed and stacked according to conventional practice and fitted 
with an additional front section to materialize the subsonic 
leading edge configuration. The above CG stacking condition 
then could be fulfilled by similarly fitting a rear section to 
restore the symmetry of the mass distribution with respect to 
the axis of maximum inertia of the profiles. This, however, 
would maximize the additional blade mass and the elongation of 
the profile chord lengths required by the compound sweep design, 
which is structurally and aerodynamically undesirable. Proper 
selection of the radius of sweep reversal thus is necessary to 
ensure minimum blade stress and aerodynamic performance penalties. 
Adjustments of the profile chord lengths can be used only to com- 
pensate for a slightly off-optimum location of the sweep re- 
versal point. Accordingly, the optimum stacking should yield 
hub stresses exceeding those of a conventional blade only by 
the contribution due to the blade mass added to Incorporate 
the subsonic leading edge configuration. 


Prom the preliminary design iterations, the meridional 
projection of the subsonic leading edge line and its sweep re- 
versal point were known with sufficient accuracy to define the 
radial distribution of the relative Mach numbers M (r) and 

the relative flow angles 3^ (r) and (r) at the leading edge 

for final design. Those data were interpolated on the stream- 
lines between stations 9> 10, 11, 12 of the R-121 flow calcula- 
tion. (For the axial station nomenclature, refer to Pig. 14 
and Appendix A.) Table 3 presents the interpolated inlet Mach 
numbers M^ , together with the selected Mach factors f and the 


corresponding Mach numbers M^ of the relative velocity component 

perpendicular to the leading edge and M^* of the relative velocl- 

1 

ties, introduced in Eqs. 1, 2 and 4. 


I 


[ 


01 


I 


38 


r 




i 


TABLE 3. Interpolated Aerodynamic Data for Final Subsonic 
Leading Edge Design 


Leading Edge 
Radius (mm) 

Relative Mach 
Nr. 

Mach Factor 
f 

1 

^ W J ^ J 


1l 

no 

.829 

1.206 

1.000 

.829 

116 

.859 

1. 170 

1.005 

. 855 

122 

. 888 

1. 140 

1.012 

. 877 

129 

.924 

1. 118 

1.033 

.894 

136 

. 959 

1. 107 

1.062 

. 903 

! 143 

.994 

1.102 

1.095 

.908 

150 

1. 028 

1. 100 

1. 131 

. 909 

160 

1 

1. 078 



1. 186 



170 

1.127 



1.240 



180 

1. 185 



1.304 



190 

1.242 



1.366 



200 

1. 302 



1.432 



210 

1. 363 



1.499 



220 

1.422 



1. 564 



230 

1.481 



1.629 



240 

1.539 



1.693 



249 

1. 588 

1. 100 

1.747 

.909 


] 


4o 


■J 








It will be seen that forward sweep starts Immediately at 
the hub by setting M* =1, i.e., by requiring that M h M, = 

.829 at the hub section. The selected values of M increase 

w 

gradually to .91 at approximately 1/3 of the span in accordance 
with the decreasing thickness and camber of the profiles, and 
then remain constant up to the tip section. 

The cylindrical coordinates in the lateral sweep angle v 
of the subsonic leading edge line are listed in the outlined 
columns of Table 4, which reproduces the input/output data of 
the computerized calculation. 

5.4 Rotor Blade Profile Definition and Stacking Procedure 

The optimum profile stacking configuration can be described 
as follows: At every blade section along the span, the CG of 

the upper blade portion projects radially on or near the axis 
of minimum inertia of that section. This means that the radial 
projections of the individual CG’s of the upper profiles must 
straddle the axis of minimum Inertia of the lower section (sub- 
sequently referred to as i-straddling) . This is achieved by 
iterative selection of the lateral sweep angle v along the span. 
During that iteration, the radial location of the point of sweep 
reversal initially selected is kept unchanged. When adequate 
1-straddling is obtained for all blade sections, the CG straddl- 
ing with respect to the axis of maximum Inertia (I-straddling) 
of the hub section is checked and the radial location of the 
point of sweep reversal modified accordingly. 

The first preliminary design investigations were carried 
out with double circular arc profiles. In the course of the 
profile stacking iterations, it appeared that using airfoil 
sections with CG’s shifting progressively backward in the lower 
span portion with forward leading edge sweep, and forward in 
the upper portion with backward sweep, i.e., a blade configura- 
tion with minimum chordwlse excursion of the profile CG’s, 
could substantially contribute to minimize bending stresses. 

A simple analytical blade thickness distribution was used 
to simplify the design iterations involving changes in section 
properties to help relieve stresses. The thickness distribution 
is written in the following parametric form 

t(x) = kx^ (c-x) (5) 


41 



TABLE 4. CYLINDRICAL COORDINATES AND LATERAL SWEEP ANGLE OF THE SUBSONIC 
ROTOR LEADING EDGE LINE. 





where c is the chord length and n, a shape parameter. By adding 
a leading and trailing edge thickness. 


where t is the LE and TE thickness factor and v = t /C the 

ni3.x^ ^ 

relative blade thickness, a practical blade thickness distri- 
bution is obtained. The abclssa for maximum thickness is given 


max 


nc 

n+1 


the factor k by 


( 6 ) 



The complete non-dlmenslonallzed formula is 


t 

c 


= T V 



(7) 


For n 1, = .5. Furthermore, the second 

constantT^^ so that the resulting profile is 
essentially a double circular arc profile for smLl thickness. 


For n > 1, (x/c) > 1/2 and the profile CG shifts 

max 

toward the trailing edge. Since the first and second deriva- 
tives of the thickness distribution are continuous, the orofile 
curvature is continuous. 


Using profiles with circular mean camber lines and n varying 
from 1 to 1.8 from the hub to the point of sweep reversal, 
and back to 1 at the tip section, a favorable blade configuration 
was obtained. However, manufacturing difficulties and the extreme 
sensitivity to tolerance and foreign object damage of thin 
profiles with n > I .5 lead to the selection of n=l, i.e., essentially 
double circular arc profiles for final rotor blade design. 
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The blade cascade geometry was defined by means of con- 
ventional procedures and criteria. Figure 19 shows representa- 
tive streamline velocity triangles, together with the correspond- 
ing relative flow deceleration rates W 2 /W 1 and static pressure 
ratios Pj/Pj, the selected cascade solidities 0 = c/s and the 
resulting D- factor values. The hub and tip cascade solidities 
are equivalent to those which would have been selected for a 
conventional design with identical rotor inlet and exit flow 
conditions. The 30 % streamline velocity triangles are represen- 
tative of the conditions at the sweep reversal section (r = 170mm) . 

The flow deviation angles 6 at rotor exit were calculated 
with Carter’s empirical formula (Ref. 13) 

6 = m(j)/V^ 

with m = 0.23 + 0.05 ^2 (circular mean camber line). For 

small camber angles (p, Eq. (8) gives unacceptably low deviation 

angles, especially in transonic cascades with shock-boundary 
layer flow interaction. A minimum deviation angle of 2° was 
arbitrarily assumed and the calculated 6- values were faired 
gradually to the minimum value toward the tip section. The 
actual profiles were defined on coaxial cylinders for the most 
part of the blade. Three profiles were defined on cones in the 
hub region to ensure a smooth evolution of the profile geometry 
toward the conical hub section. Fig. 20 shows the relative inlet 
and exit angles (I , ^ with the tangential direction and the 
deviation angles 6 used to define the cascade geometry. 

All profiles were set at a nominal incidence i = +2° with respect 
to the suction surface. The selected profile sections are in- 
dicated on Fig. 21. Table 5 lists the profile design data de- 
fining the cylindrical and conical sections unwrapped on planes 
tangent to the cylinders and cones. (While all angles are con- 
served in the development of cylindrical sections, the profile 
camber angle is reduced in the developed conical sections by 
the sector angle formed by the radii passing through the leading 
and trailing edge points.) 

The coordinates of the center of gravity of a cylindrical 
section are determined by the following simple relations: 





RELATIVE EXIT ANGLE /3o (®) DEVIATION ANGLE 



RELATIVE INLET ANGLE 















FIG 21. DEVELOPED CYLINDRICAL SWEEP REVERSAL SECTION WITH RADIAL PROJECTION 
OF PROFILE CG's AND LE AND TE LINES. 



TABLE 5. Rotor Blade Profile Data 

28 Blades (Developed Cylindrical and Conical Sections). 


Section Radius 
(mm) 

Mean Camber 
Angle ^ (°) 

Setting Angle 

if (»> 

Chord Length 
c (mm) 

Relative Thick- 
ness /2^(%) 

110/134 

28.85 

62.44* 

64. 55 

10.77 

122/140 

24.70 

56.45* 

67.91 

9.73 

136/145 

24.95 

50.93* 

74.73 

8.03 

150 

26. 10 

46.05 

85.20 

6. 

18 

160 

20. 90 

41.95 

95.20 

4.95 

170 

17.00 

38.60 

105.90 

4. 

30 

180 

13.70 

36.05 

102.30 



190 

10.90 

• 

33.95 

97.70 



200 

8.70 

32.25 

92.90 



210 

6.90 

30.75 

88.20 



220 

5.60 

29.40 

83.80 



230 

4.70 

28.25 

80.20 



240 

3.90 

27.05 

77.20 



249 

3.30 

25.85 

75.00 

4. 

00 


♦Angle between chord and tangent to the developed section circle at the trail- 
ing edge 






+ 


( 9 ) 


and 



. 5 c cos Y d sin v 
r 



+ • 5 0 sin y - d cos y 


( 10 ) 


where c is the chord length, and d is the distance of the CG to 
the profile chord in the developed section. 

Figure 22 shows the situation for a developed conical 
section. From the aerodynamic design, the geometric charac- 
teristics of the profile, especially the inlet and exit angles 
and between the tangent to the mean camber line at LE 
g g 

and TE and the circumferential direction, are known. Also known 
are the inlet and exit radii r^ and r^ and the meridional pro- 
jection Cm of the chord. Hence, from similar triangles in the 
meridional plane: 


— Further, + c 

and with m = R^ sin ip and 6R = R^ (1 - cos 

c^ = (Cm+6R)^ + m^ = cjj + 2R^R2 (1 - cos ’I') (11) 

In the developed section, the camber angle is 

^ = ^2 " ^1 “ ( 12 ) 
g g 

and the setting angle is defined by 


sin Y “ /c 


( 13 ) 


Assuming a circular mean camber line in the developed 
section, 

62 = Y + */2 

PC 


( 1^0 


^9 




FIG. 22. DEFINITION OF CONICAL BLADE SECTIONS. 
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Equations (Il)-(l4) determine the four quantities c, $ 

and Y • They must be solved by successive iterations. Assuming 

tentatively ip, equation (11) gives c, equation (12) gives 

equation (13) gives Y , while ip is Iterated until equation (l4) is 

satisfied. 

After a profile is superimposed upon the circular mean 
camber line, CG distance d is known and the coordinates of the 
center of gravity are determined as follows: 

2^22 /d \ 

^ = 1 J- + d (symmetrical profile), e= sin“^( /^jot = 90 + y-'4 

Hence, ® - 2 cos a and from triangle 0-LE-CG: 

sin =A sin a/R^g 


* 2 

Finally, r^g = ^ R^^ and the cylindrical coordinates of 
the center of gravity are 

R il; 

eog - e + -££_££ 

eg 

^cg " ^^cg ■ ^ (16) 


All CG stacking Investigations, including preliminary bending 
stress evaluations, were carried out manually. However, as will 
be discussed later, verification of stress levels was carried out 
using computer programs at BBN and AVCO Lycoming. Figure 23 
shows the final stacking of the profile CG’s radially projected 
on the conical hub section, which was investigated by NASTRAN 
analysis. The corresponding distribution of the lateral sweep 
angle v is shown on Table 4. The NASTRAN results indicated that 
the stress distribution at the hub section could be improved 
by a slight tangential shift of the first two conical sections 
in the rotation sense. Ae - shifts of -.008 for the hub and 
-.004 for the next section were effected without readjusting 
the z - coordinates of the leading edge points. Those shifts 
are indicated on Fig. 23 . Provision has been made in the i - 
straddling to generate a moment that continuously compensates the 
moment of the aerodynamic forces, (which are reflected in results 
hereafter) . 
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FIG. 23. CONICAL HUB SECTION DEVELOPED ONTO PLANE TANGENT TO 
CONE, WITH SUPERIMPOSED RADIAL PROJECTION OF PROFILE 



The optimum radial distribution of the lateral sweep angle 
V Is different for forward and backward leading edge sweep 
directions. Consequently, a discontinuity of lateral sweep may 
occur at the sections above and below the point of sweep reversal. 
This, In turn, results In a high rate of curvature of the blade 
surface. Since the blade Is defined by discrete sections, this 
appears only as a more or less pronounced concentration of the 
spanwlse curvature of the blade surface In the sweep reversal 
region. Nevertheless, this local curvature Increase generated 
prohibitive stresses near the trailing edge In several preliminarily 
generated configurations. 

This problem was compounded by the additional bending moment 
around the I-axls of the section of sweep reversal, due to the 
rearward location of the CG of the upper blade portion with back- 
ward leading edge sweep. The difficulty Increases since the sweep 
reversal was selected Initially so as to minimize that moment and 
It was gradually moved Inward from rg„ = l88 to 170 mm, still 
leaving the blade CG In a forward position with respect to the 
I-axls of the hub section. The stress concentration problem at 
the sweep reversal section was solved by means of an elaborate 
compromise of the profile stacking through that section. Involv- 
ing especially the selection of the critical lateral sweep angle 
discontinuity. For the final configuration, with rgp = 170 mm, 
this was achieved at a late design state only, the last optimiza- 
tion step, which would have required the sweep reversal point to 
be set at 160 mm radius, or the profile chord lengths to be In- 
creased In the upper blade section. With the present stacking, 
the highest stress Is 645 N/mm^ (93.5 ksl), which Is adequate 
for concept demonstration purposes. Figure 21 shows the developed 
sweep reversal section, together with the radial projections of 
the profile CG’s of the upper blade portion and the leading 
and trailing edge lines. The upper profile CG's have been stacked 
to compensate for the aerodynamic moment and to minimize the 
additional TE tensile stress resulting from the rearward CG 
position of the upper blade portion. 

Whereas the radial projection of the leading edge points 
Indicates a smooth subsonic leading edge line, the trailing 
edge line does not appear to be as smooth as desirable. For 
manufacturing the blade was defined by flat sections generated 
from the blade configuration defined In the cylindrical 
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coordinates used for the stacking investigations. Any minor ir- 
regularities of the trailing edge were smoothed out by a slight 
increasing of the chord lengths of a few local sections. All 
profile data are listed in Table 5* 

5.5 A Review of the Rotor Blade Design Iterations for Stress 

Optimi zation 

The main objective of the preliminary design effort (see 
Pig. 2^) was to define a stacking configuration that maintains 
the subsonic leading edge concept, i.e., satisfies the acoustic 
rotor design requirements with as low a blade stress level as 
possible. A target design goal of 725 N/mm^ (105 ksi) maximum 
steady state stress was sought for the design speed of 18,^50 rpm. 
For the selected titanium blade, such a stress level is considered 
adequate for the demonstration purposes of this program. 

As a first step in each Iteration , both manual and 
computerized beam-type stress computations were carried out to 
develop a feel for the iterative stacking procedure and to 
ensure numerical agreement. The standard AVCO Lycoming blade 
stress computer program which was used treats the blade as a 
twisted, rotating cantilevered beam with variable section 
properties, and takes into account the shroud and aerodynamic 
forces and the centrifugal restoring moments. All trial blade 
stacking iterations were analyzed with this program. 

Simultaneously, a quick, inexpensive and efficient finite 
element analysis was used at BBN to verify the results of blade 
iterations. The program, based on SAP, was operated in con- 
junction with a blade geometry generator which was based on the 
family of blade profile shapes , described previously by Eqs. 5-7, 
which reduce to a minimum the number of parameters required to 
specify a blade shape; namely, the leading and trailing edge 
coordinates, the section setting angle and camber, and the pro- 
file shape parameters. The program was therefore very well 
suited for Iterative design studies. The purpose of the simul- 
taneous effort was to provide further verification of the beam 
and manual analysis and to help identify stress concentration, 
which are neglected in the beam-type stress analysis program 
and in the manual calculations. These efforts were deemed 
necessary because the blade configuration differs radically 
from more conventional designs, and it was uncertain whether 
conventional design methods would be sufficiently accurate. 



MAX CALCULATED STRESS (KSI) 





A NASTRAN stress analysis program was used by AVCO Lycoming 
on design iterations which were considered particularly Important, 
and for the final stress computations verification. 

The evolution of the maximum blade stress levels as the 
blade design evolved through the series of trial designs is 
shown in Fig. 2H. The results of the first design substantiated 
the impractical stress level of a blade with simple forward lead- 
ing edge sweep. The initial sweep reversal radius (SR) was 
selected at 195*2 mm. The stacking for trial design 2 was such 
that the center of gravity of each of the 13 cylindrical blade 
sections used to define the blade projected radially down onto 
the axis of minimum inertia (1-axls) of the airfoil section im- 
mediately below. For Iteration 3, all section CG’s were pro- 
jected onto the 1-axls of the hub section. For Iteration 5, all 
section CG's above the sweep reversal section were projected onto 
the SR section 1-axis, while the stacking of Iteration 3 was kept 
for the lower blade sections. As can be seen, the resulting mis- 
alignment of the upper blade portion with respect to the hub 
section produced higher hub stresses. However, this design also 
showed the lowest stress level for the upper blade portion. 

For Iteration 6, the sweep reversal radius was lowered and 
the misalignment was corrected by introducing a discontinuity 
of the lateral sweep angle, (l.e., the angle between the sweep 
direction and the radial plane containing the relative inlet 
velocity), at the point of sweep reversal. By varying this 
parameter, a number of stacking combinations involving individual 
compromises within the upper and lower blade sections, were 
Investigated. Iteration 7 shows the best result obtained with 
this stacking concept. 

With the stress level still substantially beyond the 
preliminary design goal of 105 ksl, a detailed investigation of 
the stress pattern in design 7 was performed using the NASTRAN 
stress program. The excellent correlation which was obtained 
substantiated the beam-theory analysis method as a useful 
approach to analyze blade stacking changes. 

Subsequent Iterations were conducted with the optimum stack- 
ing concept described in Sec. 5.4. This stacking satisfies the 
condition that, at every section along the span, the CG of the 
entire blade portion above the section projects radially onto 
the i-axis of the section. As shown by Iteration 8, this re- 
duced the maximum stress level very nearly to the preliminary 
design target value. 


56 



The new stacking concept confirmed the necessity of a lateral 
sweep angle discontinuity at the point of sweep reversal to 
achieve proper stacking of the profile CG’s across that section. 
This discontinuity, however, resulted in a rapid change of the 
spanwlse curvature of the blade surface in the trailing edge 
region, which in turn results in a local stress concentration 
that was not shown by the simplified analysis. Iterations 1-8 
were conducted with double circular arc profiles (DCA). Itera- 
tions 9 and 10 used new profiles featuring rearward CG shifts 
from the hub to the section of sweep reversal, and forward CG 
shifts from that section to the tip (see previous section). In 
this way, the CG excursions from a radial line were minimized 
within the leading and trailing edge envelope and the stresses 
were reduced to the target level. 

Figures 25 and 26 show the moments about the axes of minimum 
inertia and maximum stress distributions for Design 10 as calcu- 
lated by the standard blade stress program. The influence of 
aerodynamic loads and centrifugal restoring moments are also 
shown. (Design 10 was chosen for further study since this is 
the design which first indicated stresses below the design goal.) 

A detailed Investigation of Design 10 was also performed 
with the NASTRAN program. The results showed local high stresses 
of 96 ksl at the trailing edge of the sweep reversal section and 
110 ksi at the leading edge of the hub section. By slightly 
increasing the chord length of the sweep reversal section, and 
slight re-alignment of the conical hub, these stresses were 
brought down to 84 and 96 ksi, respectively. The NASTRAN finite 
element representation of this configuration, called Design lOA, 
is shown in Pig. 27. The stress distributions of the suction 
and pressure surfaces are shown in Pig. 28. 

During the entire Iteration process, it was apparent that 
the radial location of the point of sweep reversal would have 
to be moved substantially inward from its initially assumed 
location in order to avoid a large moment about the I-axls 
of the hub section. Moving the point of sweep reversal Inboard, 
however, increases the bending moment about the I-axls of the 
sweep reversal section, thereby increasing the tensile stress 
at the trailing edge of that section. To minimize the local 
trailing edge stress concentration the radial location of the 
sweep reversal section was moved inboard cautiously. Even so, 
the blade CG remained ahead of the hub section I-axis, and 
resulted in an additional bending stress (on the order of 120 
N/mm*^) at the hub section leading edge. 
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FIG. 26. MAX STRESS DISTRIBUTION [Preliminary Design 10] 




FIG. 27. NASTRAN ANALYSIS [Preliminary Design lOA]. 




FIG. 28. EXAMPLE OF INTERIM RESULTS OF NASTRAN STRESS 
ANALYSIS MAX SHEAR CRITERION [Preliminary 
Design lOA]. 
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Prior to the selection and analysis of the final blade 
design, several intermediate designs were investigated based 
on local shifts of the CG location within the Indlficual air- 
foil profiles. (Noted as Designs 11 through 16 in Pig. 24.) 

The polynomial blade sections had been evolving toward n*l or 
a DCA profile. For manufacturing reasons, however, double cir- 
cular arc profiles were specified for the final design. This 
raised the stresses to virtually the level of iteration 9, and 
additional stacking iterations were required to achieve the 
design objective. In particular, the relative blade thickness 
was increased from 10.00 to 10.7755 at the hub section. This 
resulted in a 1055 decrease in the stress level. Additional 
reductions were achieved through a judicious balancing of .the 
profile stacking in the lower blade portion and lateral sweep 
angle discontinuity at the sweep reversal section. 

A check was performed to see if the DCA profiles allowed 
adequate flow area margins. On an average basis, the rotor 
throat passage area has a large margin to sonic throat area 
because of the comparatively high mean relative inlet Mach 
number level = 1.33 and the positive inlet incidence of 2° 

selected for optimum blading efficiency. The throat hub region 
is most susceptible to local throat choking because of the 
transonic inlet flow conditions and the higher relative blade 
thickness. Because of unknown 3-dimensional flow effects, it 
is difficult to determine local blade stream tube areas and no 
definite section throat area margins thus were specified for 
the design. A check, however, was tentatively made for the 
rotor hub section. On the two-dimensional basis of the developed 
section of Pig. 23 the ratio of throat to inlet passage width is 
1.045. At the throat location, however, the channel height has 
decreased from 139 to 136.3 mm. Assuming that all individual 
stream tube heights are reduced in the same proportion, the 
effective geometric throat/inlet area ratio thus is A mi n/Ain = 
1.045 ^ 136 . 3/139 “ 1 . 027 . With a relative inlet Mach number 
of . 825 , the sonic area ratio Ain/Ag is 1.0285, thus A« 4 „/Ae = 
1.027 1.0285 = 1 . 055 , l.e., a 5 . 5 % choke area margin. 

In the hub region, the flow has the tendency to be deflected 
inwards because of the forward leading edge sweep. On the other 
hand, the increasing density toward the tip at rotor exit com- 
bines with the essentially constant axial velocity of free-vortex 
flow to shift the streamline pattern outwards at rotor exit. 

Those compensating effects cannot be quantified at the throat 
location and the comparatively large 5.555 margin thus was judged 
adequate to account for the possibility of unfavorable three- 
dimensional effects and for the suction side boundary layer 
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growth upstream of the throat in the absence of a detached lead- 
ing edge shock. In summary, in spite of the selection of DCA 

profiles, the individual rotor section throat margins are ade- 
quate. 

5.6 Final Rotor Blade Stress Analysis 

The stress analysis for the final design iteration was 
performed using NASTRAN. The loads considered in this run 
were based on the maximum operating speed of 18,^50 rpm. In 
addition to the major contribution of the centrifugal load, 
aerodynamic gas pressure loads, the centrifugal load and the 
torsional restraints of the part span shroud were applied to 
the blade. The resulting von Mises effective stress patterns 
over the pressure and suction surfaces of the blade are shown 
in Figs. 29a and 29b. An Independent verification of these 
results was performed using the SAP program at BBN. 

The maximum stress level of 645 N/mm^ (93-5 ksi) is at the 
root near the leading edge on the suction surface. The high 
stress region of 90 ksi, however, extends only over a small 
portion of the suction surface (Fig. 29b) and so should not pose 
a problem for the planned test program. The permissible number 
of start/maximum speed/stop cycles is approximately 500, con- 
sidering a notch condition (SCF = 3*5) at the juncture of the 
blade airfoil and the base shroud. 

The tendency of the blade to untwist at the shroud location 
is small since there is only 1/2 degree difference in untwist 
between the shrouded and unshrouded NASTRAN results. The most 
significant load on the shroud, therefore, is the bending load 
due to the centllevered mass. The maximum shroud stress of 
78.7 ksi is at the blade-shroud juncture, and is conservative 
in that the large fairing radius at the juncture was not in- 
cluded in the calculation. Because of the constraining effect 
of the mid-span shroud, the untwist of the blade at the shroud 
location is negligibly small. The untwist of the tip section 
calculated from the NASTRAN results is .36°, thus Increasing 
the tip incidence from 2 to 2.4° at the design speed, a value 
well within the blade incidence design tolerance. However, 
radial growth of the shroud has not been accounted for and, if 
such growth occurs, undesirable increased tip incidence angle 
could result, due to the consequences of shroud sections "un- 
locking" . 
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The magnitude of the stresses In the fan blade airfoil are 
acceptable for an experimental program. The computed stresses 
in this design exceed AVCO Lycoming practice for titanium blades 
for longtime service operation, but fall within acceptable 
limits for the planned experimental program. 

5.7 Rotor Blade Vibration and Flutter 

The avoidance of large amplitudes of resonant vibration of 
the rotor blades over the full operating range is necessary to 
ensure the structural Integrity of the fan. The design procedure 
included an assessment of the natural frequencies of the rotor 
blade so that the forced vibration response is minimized, and the 
self-excited response is eliminated. The design goal for the 
minimization of forced vibration is ensuring that the rotor blade 
cannot resonate with the first three rotational orders of excit- 
ation due to possible inlet distortions. Although higher ex- 
citation orders will exist in the Intake, it is considered that 
these levels will be minimal in the clean Inflow expected in the 
acoustic test facilities and, thus, they will not generate sign- 
ificant resonant stress levels in the blade. The avoidance of 
self -excited blade vibration flutter is mandatory, since the 
associated stress levels usually lead to blade failure in a very 
short time. The two flutter phenomena that were considered in 
the design are subsonic positive stalled flutter at part-speed 
operation and supersonic unstalled flutter at design speed. The 
criteria for avoidance of these flutter conditions are based on 
extensive experience by the engine manufacturers and are ex- 
pressed in terms of a reduced velocity parameter: u/bw, where u = 
air velocity over the blade ( m/sec), b = blade semichord ( m), 
and 0 ) = frequency of vibration in the flutter mode (rads/sec). 

The empirical design limit values for this parameter under positive 
stalled flow are 6.7 and 2.4 for the first bending and first tor- 
sion modes, respectively. The supersonic unstalled flutter de- 
sign limit at first torsion frequency was : 


u 

bo) 



< 1.05 , where M = Mach number. 


The coefficients are calculated at 3/4 span. (Since supersonic 
unstalled flutter usually occurs in vibration modes which are 
predominantly torsional, only this mode is considered.) 
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A free-standing blade, assumed fixed at the base, was used 
in the calculation of the resonant frequencies. The natural 
frequencies for the unshrouded blade are shown in the excitation 
diagram of Fig. 30 . This design is clearly unsatisfactory since 
the natural frequency of the first bending mode has a second 
order resonance in the operating speed range. The stall flutter 
coefficients are and 1.53 for bending and torsion, respect- 

ively. These values are within the safe limits which were 
established as design criteria. The supersonic unstalled flutter 
parameter is I.I 6 and exceeds the safe upper limit. 

A partspan shroud is required to raise both the first bending 
and torsion natural frequencies and avoid forced and self-excited 
vibrations (flutter). As a physical model, the shroud was assumed 
to restrict the blade motion to a uniform translation at three re- 
presentative points. 

The design analysis was' checked by mounting two spare blades in 
a fixture which clamped at the root and partspan shroud locations. 

An acoustically coupled exciter was used to vibrate the blade so 
that the frequencies and mode shapes could be obtained. The com- 
parison between the measured and theoretical static frequencies 
shown in Pig. 31, is considered good, especially In view of the 
unusual blade shape. The "measured" frequency line in Fig. 31 
is actually the theoretical centrifugal stiffening line originat- 
ing at the measured static frequencies of the first three modes. 

Figure 31 shows the excitation diagram and calculated and 
measured frequencies for the final airfoil with the partspain shroud 
located at 645? of the span (201 mm radius). The first bending 
natural frequency has been raised so that it clears the first three 
excitation orders in the operating speed range. The fourth excit- 
ation order of the first mode, (e.g., four equally spaced front 
struts) however should be avoided. Based on the measured fre- 
quencies, the stalled flutter coefficients are 1.5 and 1.0 for bend- 
ing and torsion, respectively. The supersonic unstalled flutter 
coefficient is .75. These values meet the design criteria. The 
excitation diagram shows that the torsion and bending modes are 
not coincident in the operating speed range.. This ensures that 
the modes are decoupled. 

Strain gauges will be used during the test program to ensure 
that safe steady and vibrating stress levels are not exceeded. 

In order to locate the strain gauges appropriately, a vibratory 
stress survey was conducted using strain gauges during the static 
vibration tests: Pig. 32 shows the results of this test, normal- 

ized for each mode. The vibratory stress distributions, shown as 
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DURING STATIC VIBRATION TESTS ON BLADE S/N 17. 



the combined steady and alternating stresses In the blade, are plot- 
ted for each mode In conjunction with calculated steady stresses 
at each gauge location used. Figure 33 shows the Goodman dia- 
gram for the blade material and the vibrating stresses measured 
In each mode proportioned for the most critical location. From 
this diagram It Is seen that location *F’ Is the most critical 
location In terms of combined stress In the first mode of vibra- 
tion. Location *C* Is seen to be the most critical for the 
second and third modes of vibration. It Is therefore recom- 
mended that strain gauges at positions *C* and ’P’ are used to 
monitor the steady and vibrating stresses during the rig 
running. 

5.8 Attachment and Disk Analysis 

The fan disk stresses were computed by a Lycoming finite 
element program which evaluates the loading variation through- 
out the disk accounting for the effects of rotation, tempera- 
ture gradients and elastic-plastic conditions. 

Low cycle fatigue (LCF) life was evaluated for the signifi- 
cant regions, l.e., the disk serrations, the bolt holes and 
the disk bore, utilizing statistical minimum fatigue property 
data for Tlmkln 17-22AS material. The stress/straln ranges 
utilized In the life evaluation are the stabilized values cor- 
responding to start/stop excursions to 18,^50 rpm design 
speed. 

Stress concentration factors (SCF) were evaluated for those 
areas of the disk containing a high stress gradient, l.e., the 
serrations and bolt holes. This was accomplished by ratlolng 
the peak stresses determined by finite element analyses with 
the nominal stresses In each of the two regions. 

Nominal radial and tangential stress distributions for the 
fan disk are shown In Pig. 34, while the nominal stresses In 
the serration are shown In Fig. 35* The finite element models 
for the bolt holes and serrations are handled separately. 

Stress distributions about the disk bolt-holes are given 
In Pig. 36 from which an SCF of 2.06 was calculated, so that 
the resulting LCF life Is In excess of 100,000 ”start/stop" 
cycles based on the material S-N data of Pig. 37. It has been 
concluded that the disk bore also has a calculated life of at 
least 100,000 cycles. 
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Blade root attachment stresses and corresponding material 
properties are also summarized in Pig. 35. The bending effects 
in both the root and tenon have been included. The axial width 
of the base-shrouded dovetail root was determined by a permis- 
sible bearing stress of 420 N/mm^ (6l.O ksi). This is less than 
the compression yield limit, yet somewhat beyond the level at 
which fretting can occur under prolonged operation, but which 
should be satisfactory for a limited experimental program. 

Prom a disk serration finite element analysis, the SCP was 
calculated to be 3.33 which is consistent with values measured 
from photo-elastic analyses of similar blade root configura- 
tions. The corresponding LCP life is 24*000 "start/stop” cycles 
based on the appropriate curve of Pig. 37. These fatigue lives 
are ample for the anticipated program of testing. 



SECTION 6 

DETAILED STATOR DESIGN 

This section describes the detailed design of the fan stator 
which embodies the stator noise reduction concept described 
earlier. The stator uses vanes with varying sweepback angle to 
meet the criterion of a constant subsonic rotor wake trace speed 
along the stator vane span. The use of circumferential vane 
skew (lean) was avoided primarily to simplify the manufacturing 
problem. The stator vane number was chosen to cut off the 
radiation from residual sources due to end effects In the hub 
region of the vanes at blade passage frequency. The corresponding 
residual sources at the tip cannot be cut off because the spinning 
speed of wake disturbance pattern Is supersonic at the tip. 

To determine the proper vane sweep angle distribution, the 
rotor wakes were assumed to be convected with the mean flow. 

The spatial location of the wake centerline surfaces could then 
be computed from the mean flow properties by Integration down- 
stream from Initial points on the rotor trailing edge. Since 
the rotor wake pattern spins fixed with respect to the rotor, 

It Is possible to find leading edge lines whose shape Is such 
that their point of Intersection with the rotor wake centerlines 
travels at constant speed. Moving medium effects were taken 
Into account In the actual calculation of a vane leading edge 
shape (see Appendix C for details). The trace speed was made 
constant and subsonic relative to the local flow velocity vector 
at all points on the vane span. The stator vane sweep distribu- 
tion was designed to have an effective spanwlse trace speed cor- 
responding to a Mach number of 0.8 for the traveling load distri- 
bution. 

The fundamental acoustical analysis which underlies the 
stator design concept Is presented In Appendix C. In the re- 
mainder of this section, the methods for determination of the 
vane leading edge shape, and vane number are described, and 
the aerodynamic design considerations for the stator are re- 
viewed. 

6.1 Acoustic Aspects of Stator Design 

The major noise producing mechanism of the stator Is the 
Interaction between the stator blades and the wakes shed by 
the rotor. This Interaction causes fluctuating lift at the 
stator blades; the fluctuating lift In turn can be a potential 
source of noise. The fluctuating lift Is restricted essentially 
near the stator blade leading edges (SBLE); this fact Is made 



abundantly clear from the analytical work of Pllotas (Ref. 10). 

It is also well known (e.g., Lighthill, Ref. 15) that any fluctua- 
ting lift, whether at the leading or trailing edge, whether 
acoustically compact or not, whether in a stationary or moving 
acoustic medium, acts as a dipole source of sound. 

However, irrespective of the nature of the sources (l.e., 
whether monopole, dipole or quadrupole, etc.), there are certain 
aspects of acoustics of stationary and uniformly moving media 
which need to be considered before approaching the specific task 
of stator design and related acoustic problems. Discussion of 
these fundamental aspects is provided in Appendices C.l and C.2, 
and their application to the stator design of this fan is 
described below. 


6.1.1 Criteria for non-radiation 

Acoustic wavelengths at rotor blade passage frequency are 
small compared to the stator blade span. In this case, the 
criterion for non-radiation due to unsteady forces is that the 
trace phase velocity of the force disturbance be subsonic 
relative to the local gas flow. Skewing, or sweeping of the 
stator blade, increasing the separation between rotor and 
stator, and shaping the rotor blade are techniques which can 
be used to reduce the phase trace speed. 

Proper modification of leading edge profiles can reduce the 
phase trace speed along the leading edge and also the relative 
angle between that velocity and the local flow. Both effects 
are important as it is the trace velocity relative to the local 
gas-properties flow which must be kept subsonic. 

Each individual wake shed by a rotor blade suffers a lag 
in the circumferential direction. The net effect of this lag on 
the nature of Impingement of the wake on an unswept SBLE is 
that the wake hits the SBLE at the hub first and the Impingement 

process propagates radially outwards towards the SBLE tip with 
a spanwise varying phase or trace velocity c^(r). Sweeping 
back the SBLE enhances this phase lag effect, in the sense that 
the spanwise trace velocity of wake Impingement is reduced. A 
criterion along the lines of Eqs. (C.62) and (0.63) is used to 



guarantee that the wake trace Mach number mQ is less than mu 
everywhere along the swept back SBLE. The nature of the wake 
phase lag and the calculation of the SBLE sweep angle is de- 
scribed more fully below. The successful analysis of a rotor 
wake tracing along the leading edge of a stator requires under- 
standing of a set of transformations between stator-fixed co- 
ordinates and moving medium coordinates. The derivation of the 
trace velocity in stator-fixed coordinates, and subsequent 
Gallilean transformations to gas-fixed coordinates is given in 
Appendix E. 

6.1.2 Estimate of rotor viscous wake 

Estimates of the magnitude of the rotor viscous wake at the 
leading edge of the stator have been made. The method of estima- 
tion Involved modeling the rotor blade wakes as the wakes behind 
isolated airfoils. The method is somewhat crude, as it ignores 
the interference between wakes and the axial pressure gradient. 

The variation of angle of attack at the stator which results from 
the estimated velocity fluctuations is as much as 10 degrees 
from the mean. Experience with axial flow turbomachinery wakes 
indicates that the estimated rotor wake amplitudes at the stator 
leading edge are likely to predict higher resultant angle of 
attack fluctuations than will exist in the actual rotor wake. 

This is due to the higher rate of decay of rotor blade viscous 
wakes in turbomachines when compared to isolated viscous wakes 
in free flow (see, for example, Lakshminarayana and Raj; Ref. 16). 

6.1.3 Computation of rotor wake distortion 

Contours of constant phase for rotor wakes at different axial 
locations were computed by use of a stepwise integration of the 
phase lag of the wake relative to a point in the rotor, as a 
function of radius. Cylindrical helical flow was assumed (radial 
flow velocity was assumed to be zero). Axial and tangential 
velocities used in this calculation were provided by AVCO’s 
aerodynamic design program (Appendix A). Contours of constant 
phase calculated at several stations downstream of the rotor 
are shown in Fig. 38 (see Fig. l4 and Appendix A for Station 
Locations). Contours of constant phase versus axial location on 
cylindrical surfaces, shown in Fig. 39 and contours of constant 
phase in the axlal/radlal plane, shown in Fig. 40, were derived 
by cross-plotting from Fig. 38. 


1 


79 









PHASE vs AXIAL LOCATION ON CYLINDRICAL SURFACES 












Figure 4l illustrates. In stationary coordinates, the flow 
and blade motion geometry and the equations used in the calcula- 
tion of the constant phase contours. The two terms 0(a) and 0(w) 
in the calculations represent the angular translation of a fluid 
particle and the angular rotation of the rotor, respectively, in 
the time required for the fluid particle to flow from axial Sta- 
tion 1 to Station 2. The trace velocity, relative to local flow, 
for a number of constant-sweep-angle stators is shown in Pig. 42 
The very high trace velocities near the tips result from reduced 
wake "windup" in that region, thus illustrating the limited ef- 
fectiveness of constant angle swept stators. 


6.1.4 Mach .78 leading edge stator 

A blade leading edge sweep profile for trace speeds less 
than Mach 0.8 was developed for the final rotor and flow path 
design using an Iterative method to achieve a nearly uniform 
trace velocity. The blade has a minimum sweep angle of 25 degrees 
at, approximately 1/3 of the span from the root. Sweep at the 
root and tip are 30 and 40 degrees, respectively. Figure 43 shows 
the sweep profile as well as the trace and acoustic speeds as a 
function of radius. 

The calculations assume a rotor blade reference axis at the 
axial location 12-1/2 mm. forward of the root at Station 13. It 
also assumes a stator leading edge which is radial, when project- 
ed in the r, 0 plane, and has its root 12.5 mm. downstream of 
Station 15. 

The inflow-induced radiation from an array of such variably- 
swept stator vanes is now restricted to the tip regions when 
discontinuities occur. The limitation of such radiation depends 
upon proper choice of the number of swept leading edge stators, 
which in turn depends upon the rotor blade number, rotation speed, 
and moving medium acoustical considerations. The computation of 
vane numbers is discussed below. 
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6.2 Analysis For Determination Of Number Of Stator Blades 

section determines the appropriate minimum number 
V of stator blades so that the acoustic noise radiated from 

the stator at the rotoi* blade passage frequency f^ is minimized. 

Since the swept back SBLE derived in the above section is 
of finite extent, the end effects at the SBLE hub and tip from 
the wake/SBLE interactions remain as potential sources of noise. 
The aim here is to seek a partial oiraumferential cancellation 
of these end sources. Thus, two discrete circumferential arrays 
exist, one at the SBLE hub and the other at the SBLE tips. 

Since the circumferential phase velocity c is higher at the 
tip than at the hub, one concentrates on tfte discrete circumferen- 
tial array made up of uncancelled sources at the SBLE tips. 

Also, as discussed in Appendix C, the discussion of a discrete 
array must be limited to only one frequency Wq. Choosing u 
to correspond to the fundamental rotor harmonic , (l .e . , to ° 
the rotor blade passage frequency f ^ ) , one obtains 


0) = J2B 

0 


(17) 


where is the shaft rotation in radlan/sec ( 1 19^0 rad/sec) 
and B is the number of rotor blades (B = 28). The blade passage 
frequency f^ is given by 

^ t 8600 Hz (18) 


For the circumferential phase velocity c^ 



nrt 


(19) 
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where is the radius at the stator tip (r. : 0.2i|m :: 0.79 ft). 
The corresponding Mach number ra^ is therefore given by 





1.27 


( 20 ) 


With reference to results of Appendix C.2, (Eq. C.60), mj is 
the Mach number of the gas flow parallel to the array, and m 
is the gas Mach number normal to the array. Since the array^ 
under consideration is oriented circumferentially, the gas Mach 
number m in the circumferential direction plays the role of 
mi and the gas Mach number m in the axial direction plays the 
role of m^; the radial gas ^Mach number, normal to the duct 
walls, is to a good approximation zero. Thus, we have 


m = m ~ 0.353 

1 W 


( 21 ) 


m = m s 0.582 
r & 


( 22 ) 


Note that m,. Is directed the same way as the shaft rotation 
n or the phase Mach number m . Hence, first one would like to find 
from Eqs. C.62 and C .63 whether mQ < m^ - one of the two 
necessary conditions for no radiation to occur. Substituting 
the quoted values in Eq. C.62, one finds that 


“o'" ”u> for I < O i J . (23J 

in other words, the condition for no radiation is satisfied for 
angles a that are sufficiently remote from the axial direction. 


In order to satisfy the second condition for (partial) 
cancellation of radiation from the fundamental rotor harmonic, 
it is required that 



y 


(24) 


2ir 

d. 


? ijL 
'r 


(1-mfcos^a) 
s 

(1-m^-mf cos^a) 

C CL 


where d. Is the circumferential spacing between two adjacent SBLE 
tips. ' The spacing, d- , is related to the number of stator' blades 
by the relation ^ 


Pttt* 



(25) 


The right hand side of Eq.(24) is the radiation span (ic„ - k ) 
obtained from Eq. C.60, where 2ir/X has been ^o+ ^o- 

substituted for k , X being the Scoustlc wavelength at 
frequency fj.. “o Taking sound speed c in the gas to be about 
365 m/s (1200 ft/sec) , the wavelength at blade passage frequency is 


A - ^ 5 0.14 ft ' 0.043 m 

Thus, substituting Eq. 25 in Eq. 24 , the velocity is 

(1-m^cos^a) 

a 

(1-m^-mfcos^a) 

C a 



( 26 ) 


(27) 


Since the first necessary condition (Eq. 23 ) is satisfied only 
for a restricted range of angles a, it would not pay to find 
the maximum possible value of V for arbitrary a. Instead, 

Eq. 27 is evaluated for a = tt/ 4 (as a goes from tt/ 2 to ^74 
to 0, V evaluated from Eq. 27 increases), and the result is 


V > 92 


( 28 ) 


One can now examine the application of traditional analyses 
(e.g., Tyler and Sofrln (Ref. 11)) of noise generated by rotor- 
stator interaction, the anlysls that is used primarily for low- 
speed compressors H.e., analysis is based on stationary medium 
acoustics) that Involve subsonic circumferential phase speeds 
(i.e., fir. < c) and are acoustically compact (i.e., d. < X /2). 

Xf ^ r 
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An arbitrary component (say, the predominant component 
of wake velocity deficit pattern that generates fluctuating lift 
at SBLE) a(x,r,0,t) of rotor- generated flow field near the stator 
may be decomposed into circumferential harmonics as follows 


+00 


a(x,r,e,t) = 23 


n=-oo 


Aj,(x,r) 


i[ nB(e-Qt)] 


K 


(29) 


where x and r are the axial and radial locations (and for our 
case of interest denote the locations of SBLE tips) and 6 is the 
circumferential angle. 


The noise sources (in particular, the fluctuating lift I 
generated at SBLE tips) at the stator due to the nth rotor 
harmonic may then be viewed as composed of a sum of stator/rotor 
harmonics mn. A typical interaction harmonic L (s,r,6,t) may 
be written as 


where 


L (x,r,e,t) = 
mn 


L (x,r) 
mn 


i(m6-nBfit ) 
e 


(30) 


I 


m = nB + kV 


(31) 


and where k can assume arbitrary Integral values (positive, 
negative or zero). 

The circumferential phase velocity c Cr) associated 
with Eqs. 30 and 31 can be written as o mn 


c (r) = 

^o^^^rnn nB+kV » 


(32) tt 
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Similarly, from Eq. 19 , the circumferential phase velocity for 
all the rotor harmonics n is Qr and assumes the value J2r. at 
the SBLE tips. This same value is recovered for the ^ 
interaction modes from Eq. 32, for the stator fundamental 
mode, i.e., for the case k = 0. 


Figure depicts the situation in terms of these rotor- 
stator interaction harmonics. The harmonics lie at the inter- 
sections of vertical lines passing through the nB axis for 
“ 0, ±1, ±2 ... and horizontal lines passing through the m 
axis for k = 0, ±1, ±2 ... The rotor fundamental tone 
occurring at the blade passage frequency fj, Uee Eqs. 17 and 
18 .) corresponds to vertical straight lines passing through 
“ -1 (i.e., nB = ±28). Similarly, nth rotor harmonic corres- 
ponds to frequency nfj,. The fact that attention was turned 
to the stator fundamental harmonic at frequency (see Eqs. 19 
and 20 ) means that the k = 0 stator mode was examined at 

Prom Eq. 23 , oue finds, that this stator fundamental harmonic 
barely escapes radiation. From Eq. 32 , it can be seen that' the 
same situation applies to stator fundamental harmonics (i.e., k = 0 
modes) for all rotor harmonics (i.e., arbitrary n). Thus, the 
straight line in Fig. 44 passing through these k = 0 modes 
separates the radiating and non-radiating harmonics. 


The criterion of Eq. (24) was applied to prevent the next 
candidate stator harmonics (k = ±1 modes for n = ±1) from 
radiating at the blade passage frequency f {only). The 
straight line joining these k = -1, +1 mod?s thus also separates 
the radiating and non-radiating harmonics. The flow-induced 
assymmetry in radiation span along wavemmiber is reflected in 
Fig. 44 by assymmetry of radiating and non-radiating harmonics 
around m and nB axes. Incidentally, stator harmonics lying 
in the upper right and lower left quadrants of the m, nB 
plane possess circumferential velocities that are oriented in 
the same direction as shaft rotation fi, and the harmonics lying 
in the upper left and lower right quadrants possess velocities 
that are oriented in the direction opposed to Ji. 


Finally, note that the relatively high number V of 
stator blades indicated by Eq. 28 may cause design problems 
of aerodynamic nature. For example, relatively high solidity 
at the hub, particularly for the scale model fan, is unacep- 
ceptable. Therefore, a compromise number of 59 was selected 
for V. Such a choice ensures, circumferential cancellation at . 
the SBLE hub, but not at the tip. In other words, with reference 
to Fig. 44, k = ±1 modes would radiate from the SBLE tips. 


91 





6.3 Stator Aerodynamic Design Considerations 

The stator is characterized by a backward leading edge sweep 
varying from 25 to in the meridional plane. The axial spac- 
ing between rotor trailing edge (TE) and stator leading edge (LE) 
varies approximately from two to three rotor hub chords along the 
span. 


A minimum number of 59 blades was specified by acoustic 
considerations. This, in conjunction with a tentatively selected 
hub cascade solidity cr. , = 2, resulted in a chord length of 
approximately 30mm. 

The meridional contour of the stator was shown in Fig. 14. 
Radial station 17 crosses the leading edge, 19 the trailing edge, 
and 18 crosses both the leading and trailing edges. Radial 
equilibrium along those stations is markedly influenced by the 
varying degree of stator turning, resulting in peculiar tangen- 
tial velocity distributions that have been input in R-121, to- 
gether with the corresponding total pressure loss distributions. 

The flow conditions from rotor exit station 13 to stator inlet 
stations 16, 17, 18, and 19, have been calculated according to 
constant rotor exit momentum V, »r specified along the stream- 
lines. ^ 

The meridional flow pattern (Fig. 14) shows the radial 
streamline shifts induced by the swept back stator configuration, 
especially in the hub region, where is large and has a strong 

effect on radial equilibrium. In the axlsymmetrlc flow case, the 
streamlines approaching the leading edge are deflected inboard. 
Looking at the lower portion of station 17, the flow at the hub 
section has already undergone the major part of its turning, and 
V^*r thus increases markedly from the hub to the 30 % streamline 

on that station. This substantial departure from free-vortex 
flows generates an increase of the axial velocity component to- 
ward the hub and a corresponding increase of the mass flow den- 
sity PV^, in turn resulting in inboard streamline shifts between 

leading edge and station 17. This characteristic pattern is 
found along the entire span, but disappears gradually toward the 
tip section because of the decreasing value of the V^/r-term. 

The meridional streamline curvature term V^/R has a strong 

me ^ 

effect in this flow field region. The determination of R how- 

e 

ever is very approximative even with the spline-on spline procedure 
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used in R-121, so that the Interpolated values of the flow condi- 
tions at the stator leading edge cannot be expected to be smooth. 
Pig. 45 shows the radial distribution of the inlet angles 

determined by V^^-interpolation and constant V^t along the stream- 
lines, and the smooth distribution assumed for blading design. 

The maximum smoothing error does not exceed 1-1.5°, which is well 
within the accuracy that can be expected from the axisymmetric 
analysis . 

Figure 45 also shows the stator exit flow deviation angles 6 
calculated with Carter’s formula [Equation ( 8 ), circular mean 
camber lines]. The blade sections are stacked with the leading 
edge in a meridional plane. All profiles were set at 0° nominal 
Incidence. Table 6 lists the profile data defining plane sec- 
tions perpendicular to the radius in the leading edge plane. 




FIG. 45. STATOR INLET FLOW AND DEVIATION ANGLES. 
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STATOR DEVIATION ANGLE 5 (”) 



TABLE 6. STATOR BLADE DATA 



1 


I 


1 
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SECTION 7 


COMMENTS ON RESIDUAL NOISE SOURCES AND NOISE LEVELS OF THE 
SWEPT ROTOR AND STATOR FAN STAGE 

The object of the Low Source Noise Pan Program is to 
design fan components for minimal noise generation. This 
has been done by first using physical models for each of the 
component noise mechanisms, and calculating the appropriate 
parameters from the particular baseline fan design, then 
modifying the component geometry to minimize noise generation. 

All sources of noise cannot be eliminated, and indeed all sources 
have not been attacked in this study. 

7.1 Residual Sources for a Fan outside the Laboratory Environment. 

As has been previously discussed in detail, the compound 
sweep required on the rotor blades for structural reasons will 
lead to a conical shock at the sweep reversal point. However, in 
some future fan designs, the location of the sweep reversal point 
^t a radius less than that at which the critical relative Mach 
number occurs will eliminate the source of noise. Rotor 
discrete frequency mechanisms which cannot be eliminated Include 
the so-called Gutln noise sources associated with steady loads 
and thickness noise. However, the non-radlal blading may cause 
these mechanisms to excite high order duct modes and thus re- 
duce the radiation to the far field. Rotor broadband mech- 
anisms are relatively poorly understood quantitatively (in the 
absence of Inflow turbulence), and thus are difficult to attack at 
the source. Shock/turbulence Interaction in the channel may cause 
some forward radiated noise, and quite likely causes aft-radiated 
broadband noise. 

Stator noise mechanisms are much better understood and can 
be attacked with much more confidence than some rotor mechanisms. 
The uncancelled tip radiation (calculated in Appendix C) is the 
only discrete-frequency mechanism Inherently associated with 
the subsonic trace speed swept stator concept, assuming that 
the rotor wake field can be accurately specified. Stator 
broadband mechanisms not attacked by the swept leading edge in- 
clude vortex shedding and flow separation at the trailing edge. 

Other broadband noise from an installed fan comes from the 
exhaust jet and duct boundary layer turbulence interaction with 
the lip of the fan duct . 
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7.2 Prediction of Noise Levels and Noise Reduction of the 
Swept Rotor and Stator Fan. 


Despite intensive research efforts in the past twenty years 
which have. led to a good understanding of noise mechanisms and 
scaling laws, the ability to predict fan noise for an arbitrary 
design on a component-by-component basis is quite limited. For 
conventional fans, useful semi-empirical correlations of data 
have been made using scaling laws which are based on assumed 
mechanisms. Thus, for conventional fans, one can predict 
within a few dB the expected sound power and directivity. 
However, the applicability of those correlations to a fan of 
unconventional component design is doubtful. 

For the subject fan design, the prediction of residual 
noise from the rotor requires the knowledge of the strength 
of the conical shock upstream of the rotor, which is not 
presently known due to the cessation of activity on the 3-D 
compressible flow program. The stator discrete noise has been 
calculated directly for basic principles and is presented in 
Appendix C. 

However, the main noise source of Interest, rotor multiple 
pure tones cannot be reliably estimated without detailed infor- 
mation on shock structure and duct propagation characteristics. 

In the interest of providing an estimate of the benefits of 
eliminating MPT noise, a computer program published by Burdsall 
et al. y (Ref. 20) was exercised (see Appendix D for details). 

The results summarized below for a full scale (a 40,000 lb 
thrust) counterpart of the 20 inch fan built in this program, 
show that elimination of the shock-generated MPT's reduces the 
overall and perceived noise levels by 4-6 dB, and reduces the 
tonal content in the 1/3 octave band containing the blade passage 
frequency by about 10 dB'. 



TABLE 7 . ORDER-OF-MAGNITUDE EMPIRICAL ESTIMATE OF NOISE LEVELS 
FROM FULL SCALE SINGLE STAGE FAN. 


Spectrum 

Component 

Overall ,, 

PWL (dB re 10 "^^w) 

OASPL*(@ 150 ' ) 
dB(re 2 xlO"®N/m^) 

PNL*(@ 150 ;' ) 

M.P.T. (conven- 
tional blades) 

152-154 

104-106 

117-119 

B.P. Tone 

143 

95 

108 

Broadband Mech- 
anisms 

150 

103 

114 

TOTAL 




with MPT’s 

153.5-155.2 

105.3-107.2 

118-120 

without MPT’s 

150 

103 

114 


*Valid in the forward-radiated direction at 
40-80° from fan axis (± 3 dB); to scale to 
subtract 20 log r/150» where r is distance 


azimuths from 
greater distances, 
in feet . 
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SECTION 8 

MECHANICAL DESIGN ASPECTS AND FACILITY INTEGRATION 


The fan rig Is built to conventional standards and Is de- 
signed to Interface with the NASA W-2 and W-8 test facilities. 

The W-2, acoustic facility Is arranged for the measurement of 
forward and rearward radiated noise j thus, the rig casings 
have flanges at both ends which mate with the facility mounting 
flanges. The manner In which this Is accomplished Is shown 
schematically In Pig. 12. In the reverse flow mode, for back- 
ward radiated noise, an additional flow path adaptor supplied by 
NASA and not shown In Fig. 12 Is fitted to the fan outlet flanges. 
In the W-8 facility the fan Is mounted on Its rear flange with 
the flow entering from the bellmouth. All detailed performance 
measurements of the fan will be made In the W-8 facility. 

A flow path adaptor fits over the facility bearing housing 
to control the fan outlet flow and Into this adaptor Is fastened 
the Inner shroud of the stator vane assembly. The outer shroud 
of the stator vane assembly Is located In the fan casing and pro- 
vision Is made for axial adjustment of the stator by relocating 
the spacers at the Inner and outer shrouds. The fan outer casing 
Is split In the vertical plane for assembly purposes. The section 
of the outer casing In the area of the blade tips Is relieved and 
an abradable shroud lining Is Installed to prevent blade tip 
damage In case of tip rubs. Figure 46 shows the engineering 
cross-section of the fan which details all the major components. 
Strain gauges will be applied to the rotor blades, the wires 
being led down the front and rear faces of the disc . In the 
W-2 facility the sllprlng Is Installed at the driven end of the 
rig shaft and, thus, the strain gauge wiring will pass down the 
length of the hollow shaft . When the fan Is running In the W-8 
aerodynamic facility, the strain gauge wiring will be led forward 
through the drlveshaft adaptor which Is fitted In place of the 
spinner support cone. A static fairing Is Installed over the 
sllprlng to provide a smooth flow profile Into the fan. In place 
of the spinner . 




SECTION 9 


CONCLUDING REMARKS 


A research program was undertaken to try to demonstrate that 
source noise reduction concepts which are based upon full and 
rigorous application of fundamental aeroacoustlc principles 
can be implemented on turbofans in the currently-operating range 
of tip speeds and pressure ratios, utilizing the existing design 
and manufacturing capabilities of the aircraft engine industry, 
without serious compromise of the noise reduction concept. 

The subsonic leading edge rotor blade concept has significant 
potential as a practical solution to rotor-generated noise due 
to its Inherent lack of sensitivity to off -design-point operating 
conditions, and the large family of detailed edge and generating 
surface contours available for fans in various speed ranges. The 
aerodynamic behavior of subsonic leading edge rotors in supersonic 
absolute inflow velocities is largely unknown at this point in 
time. However, it is believed that the characterizations of such 
flow fields, to the extent necessary in developing actual engines 
using the subsonic leading edge rotor principle, would at this 
time require considerably less effort than has been expended 
historically in understanding aerodynamic behavior of conventional 
rotors. 

The subsonic trace speed stator vane concept can be imple- 
mented through application of moving medium acoustic principles 
and a knowledge of the details of the rotor wake field, the lack 
of the latter being a current limitation. However, the subsonic 
trace speed concept can, in principle, be successfully Implemented 
by use of conservative assumptions about the rotor wake field. 



APPENDIX A 


COMPUTER LISTING OF AEROTHERMODYNAMIC PARAMETERS 
FOR FINAL ROTOR, STATOR & FLOW PATH DESIGN 



APPENDIX A: DETAILED AEROTHERMODYNAMIC DATA 


This appendix contains a computer listing of the aerothermo- 
dynamlc data for the final design of the fan. 

The first four pages, A-6 to A-9 are input data to AVCO 
Lycoming Program R121 at various axial stations. All units 
in the SI system and headings on the columns are self-explanatory. 
The three parameters in the left hand column are: 

TOT PRESS = Total Pressure Ratio 

TOT TEMP = Total Temperature Ratio 

VU = Absolute Tangential Velocity Component of 

the Air (m/sec) 

The remaining pages are detailed output at the various 
axial stations, the non-obvious terms of which are defined 
below. . 



Coded Term 

Meaning 

Units 

A STATIC 

ambient sound speed 

m/sec 

A TOTAL 

sound speed based on total temperature 


ALPHA BAR 

sln“' (Vm/V)= angle of flow made 
by V in tangential direction measured 
on a cone 

degrees 

ALPHA 

sin ^ (Vx/V) = angle made by projec- 
tion of absolute velocity vector 
(on a cylinder) 

degrees 

BETA 

angle that the relative velocity 
vector makes with a cylinder 

degrees 

V 

absolute velocity 

m/s 

VM 

meridional component of V 

m/s 

VR 

radial component of V 

m/s 

S-VALUE 

radial length measured along a 
station cut (origin at hub) 

m 

% SPAN 

percent radial distance compared to 
full span measured from hub 

— 

VX 

axial component of absolute velocity 

m/s 

VU 

tangential component of absolute 
velocity 

m/s 

W 

relative velocity 

m/ s 

WU 

tangential component of relative 
velocity 

m/s 

MV 

Mach number of absolute velocity V 

— 

MVX 

Mach number of VX 

— 

MVM 

Mach number of VM 




A- 3 



Coded T-erm 

Meaning 

Units 

R-ADC 

streamline radius of curvature in 
meridional plane 

m 

RHO 

fluid density 

Kg/m® 

ROTOR EPP 

polytropic rotor efficiency 


S-VALUE 

radial length measured along a 
station cut (origin at hub) 

m 

STAT PRESS 

static pressure 

bars 

STAT TEMP 

static temperature 

° Kelvin 

TO/TO 1 

( TO/TO )T J 

Similar definitions for the stagnation 
temperatures . 


TOT PRESS 

total pressure ratio 


TOT TEMP 

total temperature 

° Kelvin 

U 

rotational speed 

m/s 

V 

absolute velocity 

m/s 

VM 

meridional component of V 

m/s 

VR 

radial component of V 

m/s 

VU 

tangential component of absolute 
velocity 

m/s 

VX 

axial component of absolute velocity 

m/s 

W 

relative velocity 

m/s 

WU 

tangential component of relative 
velocity 

m/s 

X-VALUE 

axial location of station re: origin 
(station 4) 

m 



Coded Term 


Meaning 


% AREA 


% of annulus area taken up by a stream 
tube from the preceding area to that 
where % AREA is indicated 
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APPENDIX B: GEOMETRIC CONSIDERATIONS FOR SUBSONIC LEADING 

EDGES ON TRANSONIC ROTOR BLADES 


We first note that a simple leading edge configuration is 
obtained by sweeping each leading edge element dl in the plane 
formed by the local relative velocity Wj and the radius (W^-r plane). 
This plane Intersects the Mach cone along two generatrices that 
form the Mach cone angle y with W. Any other plane through the apex 
cuts the cone along generatrices forming a smaller angle y" with 
W in the W-r projection. Since the radial projection of dl is 
essentially proportional to sin y”, it follows that the simple 
case defined above yields the shortest possible swept blade 
length for a given annulus height and a given relative velocity 
distribution W(r). For structural reasons, however, the leading 
edge must be swept aside from the W^-r plane. 

The general situation is shown in Pig. B.l (Refer also to 
Fig. 17). 

la shows the projection on a plane perpendicular to the radius 
passing through leading edge point P. The velocity triangle 
is projected in that plane for visualization convenience. 

lb shows the projection on the W-r plane, with the Mach lines 
forming the Mach angle y with W. In general, W forms an angle e 

” 1 

with plane la. 

Ic shows the projection on a plane perpendicular to W, Intersecting 
the Mach cone along circle c. 

Id shows the projection on a meridional plane. 
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On la, the leading edge element dl has the projection dl' and 
the radial planes passing through dl and W form the angle v, which 
is a design parameter to be selected so as to minimize the blade 
bending stresses. The resulting lateral sweep component dn 
appears also on projection Ic and causes the Mach angle y 
between dl and W to project into the W-r plane lb with a smaller 
aperture y". 


Prom la 



dx^ = ±dl' 

cos ( 3 

+v) 

dz = ±dl' 

sin ( 3 

+v) 


In the above and the following relations, the top signs denote 
backwards, the bottom signs forward sweep. 


Since 


rti « - dm _ dp 

cos V “ tan(y" * cos v 

1 


6 + v) 


la. • 

L J. 


tan (y”±e„ ) 


cos V 


dp 


Prom lb and Ic 


tan y” » ds cos (p ^ tanV * cos w 
dX 


(b.i) 


(B.2) 
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and since 


sin 0 ) = 


dn 

ds 


dn 

dm 

dl” 

COS (A) 

dm 

dl" 

ds 

cos 0) 


= tan V • cos (y"± e ) ^ 

siny" 

Therefore, 


COSb) 



tan^v 

sin^y” 


1 


cos^ (y" ±e ) 


which is introduced in Eqn. B.2, yielding 


(B.3) 


tan^y" = 


tan^y 


1 + 


tanlv 


sin^y" 


cos^ ( y" ±e..) 


(B.4) 


Developing cos ( y” -e„), Eqn. E4 is reduced to a quadratic 

W 


tan 


equation for tan y"* The solution is 

4 


±sin cos e^tan^v +ltan*y (l+sin^e^an*v)-cos* ^ 


y" = 


(Only the (+) sign is valid in front of the radical, since 
tan y” must tend toward tan y when v tends toward 0). 


We define the blade profiles and their stacking in cylin- 
drical coordinates. The angular abclssa of leading edge point 
P(r) then is 



1 1 


Eqs. (B.6) and (B.7), together with (B.5) and (B.6), determine the 
cylindrical coordinates of the profile leading edges, for a blade 
with sonic leading edge. 

With the section profile data, the stacking of the centers of 
gravity is determined, and the blade bending stresses can be 
calculated. However, it is advisable to iterate the leading 
edge coordinates until a favorable alignment of the CG’s is 
achieved, prior to the calculation of stresses. 



By optimum selection of the lateral sweep and of the 
radial location of the point of sweep reversal, it is expected 
that the additional stresses affecting the subsonic leading 
edge blade will be reduced to: 

(a) Additional centrifugal stresses from the added 
blade mass necessary to materialize the subsonic 
leading edge configuration. 

(b) Bending stresses from moments without substantial 
component in the direction of the axes of minimum 
inertia. 
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FUNDAMENTAL ACOUSTICAL ASPECTS OF STATOR DESIGN 


C.l Continuous and Discrete Line Sources in a Stationary 
Acoustic Medium 

Continuous Line Source 

Consider a line monopole source of the type 

l(koX-O)^t) 

q(x,t) = Qq e , (C,l) 

where is the source strength per unit length. The line 
source of Eq. C.l represents one wave traveling along the 
X-axis (see Pig. C.l) with a velocity Cq given by 

One is interested (only) in the far field pressure p(x,y,z,t) 
radiated by the line source. Define 

r = (y^+z^)^ (C.3) 


Consider the case (referred to as Case No. 1) where the source 
velocity c^ is supersonic, i.e., 

|Cq| > c , (C.4) 

or equivalently, 

l^ol < “ l%l/° • (C.5) 

0 

Here c is the sound speed for the medium and is the acoustic 
wavenumber at frequency ojq. For this case, the°far field pres- 
sure p(x,r,t) is non-zero; in other words the line source can 
radiate acoustic power. More specifically. 




p(x,r,t) = constant 


X 


e 


9 


(C.6) 




(k^r)’ 


where the radial wavenumber is given by 

2 \h 


*'r " <*'1 - 

o 


(C.7) 


Since k^^ < k| (see Eq. C.5), k„ is real and the sound is 

° propagated radially outwards from the x-axis. 

Now consider the alternate case (Case No. 2) where the 
source velocity Cq Is subsonic, i.e.. 


or equivalently. 

O 

o 

A 

o 

(C.8) 


l>'ol ’ *=3 • 
o 

(C.9) 

For this case, the far field pressure p(x,r,t) is 
words the line source cannot deliver any acoustic 
specifically. 

zero; in other 
power . More 


p(x,r,t,) = 0 

(C.IO) 

This happens because 

the radial wavenumber k is 

r 

Imaginary, 


■'r = "1 <>'o - • 

(C.ll) 


and the near field pressure decays exponentially in the radial 
direction. 

Let us reconsider the above results in terms of the spatial 
Fourier transform q(kj,t) of Eq. (C.l). First, define the 
general Fourier transforms . 


C-4 


dx 


(C.12) 


q(kj,t) = ^ j q(x,t) 


-IkjX 


q(x,t) = 


ik X 

q(k ,t) e ^ dk 


(C.13) 


Unless stated otherwise, the limits of integration are always 
to be taken from -® to +». For later use, the temporal Fourier 
transforms shall also be needed, defined as follows. 


= 27t 


io)t 

q(k^,t) e dt 




-l(k x-oit) 

q(x,t) e ^ dxdt 


q(x,t) = 


l(k x-(ot) 

q(k , 0 )) e * 

1 


dkjdw 


iu)t 

q(x,w) e doj 


(C.llO 

(C.15) 

(C.16) 

(C.17) 


Substituting q(x,t) of Eq. (C.l) into Eq. (C.12), 

-lo) t 

q(k^,t) = Q^e ° 6(k^-k^), (C.l8) 


where 6 is the Dirac delta function. Figure C.2 Illustrates 
q(kj,t) for Case No. 1 (radiation) and for Case No. 2 (no 
radiation). The "radiation span" along the wavenumber kj is 
centered around the wavenumber kx = 0 and ranges from -k« to 
+kg^ . This radiation span is shown in Pig. C.2 as a ° 
o shaded strip. 

Let us reformulate the condition of no radiation in terms 
of velocities and Mach numbers. The two extremes -k. and +k 
of the radiation span correspond respectively to the^o 
lowest £ind the highest velocities that the source wave can have 
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q(k,,t) THE SPATIAL FOURIER 
TRANSFORM OFq(x.t) 


NO RADIATION 
FIG. C.2. 


CASE I 




CASE 2 


-k 


Go 


*^ao 


RADIATION 


NO RADIATION 


CASES OF RADIATION (No. 1) AND NO RADIATION (No. ^) 
ILLUSTRATED IN TERMS OF THE WAVENUMBER k . 
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for no radiation to occur (Case 2). These extreme velocities 
and Mach numbers are, 


c 


Z 



m 


I 



c 


-1 


c 


u 


m 


u 



= +c 


+1 


(C.19) 


(C.20) 


The source wave Mach number Is defined as : 



(C.21) 


Thus, the condition of no radiation (Case 2) becomes 


m, 


< m_ < m 


u 


(C.22) 


Next, consider a spatially frozen but otherwise arbitrary 
pattern q(x,t) traveling, as before, with fixed velocity Cq In 
the x-dlrection. Thus, 

q(x,t) « Q(x-c^t) . (C.23) 

In contrast to Eq. (C.l), for which there was one wavenumber ICq, 
one frequency Wq* and one velocity Cq, for Eq. (C.23) there is 
a range of wavenumbers k^, a corresponding range of frequencies 
CO, and one velocity Cq. Using Eq. (C.15), the Fourier transform 
of q(x,t) of Eq. (C.23) Is 
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q(kj ,w) 


1 


dxdt (C.24) 


(2tt) ^ 


-l(k x-cot) 
Q(x-c^t) e * 


2n 


^ -Ik c t+icot 
Q(k) e ' ° dt 


(C.25) 


= Q(k) 6(o)-k c ) , 
1 o 


(C.26) 


where 


Q(k) 


2tt 


-Ik X 

Q(x) e ^ dx 


(C.27) 


Figure C.3 shows the straight lines along which q(kj ,03) of 
Eq. (C.26) is non-zero. Analogous to Fig. C.2, straight lines 
corresponding to Case 1 (radiation) and Case 2 (no radiation) 
are Illustrated. Notice that as frequency oi Increases, the 
radiation span 2k^ over wavenumber k^ also increases linearly. 
But, as long as the source velocity magnitude |cq| is subsonic, 
there is no radiation to the far field. Incidentially , the 
upper right and the lower left quadrants of u-k^ plane correspond 
to positive phase velocities, (i.e., velocities along increasing 
x), whereas the upper left and the lower right quadrants of w-k 
plane correspond to negative phase velocities. ^ 

This completes the discussion of a continuous frozen pattern 
of line sources in a stationary acoustic medium. Consideration 
of the fact that the acoustic medium is moving uniformly, will 
simply alter the radiation span along kj , as will be discussed 
in Sec. C.2. However, a frozen convecting pattern along x^, 
will or will not radiate, by exactly analogous rules as developed 
here, i.e., in terms of the convection or phase velocity Cq of 
the pattern. 

Finally, note that for a continuous convecting line source 
of finite length, even if the convection velocity c^ is subsonic, 
there will be inevitable radiation from the two ends of the line 
source. For low enough frequencies, the two ends may be less 
than half an acoustic wavelength apart, in which case there may 
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FI6i C.3. CASES OF RADIATION AND NO RADIATION FOR A SPATIALLY 
FROZEN ARBITRARY PATTERN* ILLUSTRATED IN THE k .w 
PLANE. ^ 
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N 


be partial cancellation from the two end sources. For higher 
frequencies, the two end sources will radiate independently. 

This last remark is discussed more fully below when discrete line 9 

sources are considered in a stationary acoustic medium. However, 

the SBLE is regarded as a continuous line array, and a typical 

rotor wake impinging on it has a local convection velocity Cq 

along the span of the SBLE, Thus the preceding discussion of 

continuous line sources, or rather its related extension in Sec. 

C.2, where account is taken also of the moving-medium acoustics, 

is applied to determine the SBLE sweep; the criterion that is 

applied is in terms of the spanwlse local velocity of the rotor 

wake along the SBLE. If 


Discrete Line Source 

Now consider an array of equlspaced coherent monopoles, 
spaced d apart (see Fig. C.4), where: 

Xj = dj , J = 0, ±1, ±2 ... (C.28) 

In analogy with Eq. (C.l), consider one wavenumber ko, 
one frequency Uq and' the corresponding phase velocity c^. 

Thus, the source number J has the strength q(x.,t) given by 

J 

i(k X .-0) t ) 

q(x, ,t) = Qjj e ° J ° 6(x-dj) (0.29) 

r 

The entire source strength can be written as 

+® l(k x-o) t) 

q(x,t) = Qq e ° ° 6(x-dJ) , (C.30) 

J=_oo 


and the phase velocity c^, as before, is given by 

c = w /k . (C.31) 

o 0 0 


Eq. (C.12) is used to find the spatial Fourier transform of 
q(x,t) of Eq. (C.30) , 
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[ 




Q. 


q(k,,t) = 


27T 


Q 


+ 00 


-ik^x Kk^x-w^t) £ 


Q -Iw t 

° ° E 


1 0 


j=-00 

6(x-dj) dx 


(C.32) 

)dj 

(C.33) 


j=_o 


-iw t +«> 

° e ° E 

^ O 1 


2fnn> 


(C.34) 


I 


m=-“ 


q(kj,t) thus consists of an Infinite string of Dirac delta 
functions equispaced along the wavenumber k^, the spacing 
between two adjacent delta functions being 2ir/d. It is 
only the "fundamental mode" or harmonic at k^ “ m = 0 

in Eq. C.3^) that corresponds to the trace velocity c„ of Eq. 
(C.31). The remaining infinite harmonics correspond to infinite 
other velocities. The rule of radiation (Case 1) or no radia- 
tion (Case 2) is exactly the same as the one developed for the 
continuous array and illustrated in Fig. (C.2). If the funda- 
mental mode or any harmonlc(s) lie within the radiation span 
(“k^ 3 k^ ) , radiation will occur from the fundamental mode or 
° ° from the harmonic (s) lying within the radiation span. 


However, a more interesting and new feature of the discrete 
array is the classification based on a different criterion. 

That classification is as follows: 


Case A,'- 




or d 



(C.35) 


Case B, 




or 



(C.36) 


For Case A, the spacing 2u/d along wavenumber k between harmonics 
is greater than the radiation span 2kg^ , since = 2 tt/X , where X‘ 
is the acoustic wavelength at frequency w » ° 

the same condition is expressed by'the statement that the array 
spacing d is smaller than half the acoustic wavelength. For 


D 


II 


a 
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Case B, the opposite situations occur in the wavenumber and 
spatial descriptions. 

The importance of these two cases is depicted in the next 
few figures. Figure C.5 describes Case A1 (the numbers 1 and 2 
denote the older classification, 1 corresponds to radiation 
occurring, 2 corresponds to radiation not occurring). The radia- 
tion occurs from the fundamental mode at kj = k^, but since 
2ir/d > 2k^ , no other harmonic can radiate. Figure C.6 also 
describes ® Case Al, however, this time a harmonic, and only one 
harmonic radiates. Figure C.7 shows the Case A2, a situation 
one would hope to achieve. The fundamental mode at k^ = k^ lies 
just outside the radiation span, and no harmonic lies within 
the radiation span, hence no radiation occurs. Note that for 
this desired situation, the constraint of Eq. (C.8) (or equi- 
valently of Eq. C.22) as well as the constraint of Eq. (C.35) 
applies . 

Finally, Fig. C.8 shows Case Bl. There is no Case B2. 
Radiation must occur through same mode(s), whether the phase 
velocity Cq is subsonic or supersonic. In other words, arranging 
for Case B, i.e., having array spacing d greater than half the 
wavelength, is basically a poor design. 

Note that in contrast to the continuous line array for 
which the discussion related to Fig. C.2 for frequency could 
be generalized to discussion related to Fig. C.3 for all fre- 
quencies, the discussion of the discrete array presented above 
cannot be similarly generalized to all frequencies. This is 
because the array spacing d is in general fixed, whereas the 
radiation span 2k^ increases linearly with increasing (Dq* 

Thus, Case A for o frequency Wq is bound to merge into Case B 
at some higher frequency. 

Aside from the re-deflnitlon in Sec. C.2 of the radiation 
span in wavenumber kj, induced by consideration of moving -medium 
acoustics, the above discussion of a discrete array is applied 
to determine the number of stator blades, the spacing d correspond- 
ing to the circumferential spacing between two adjacent stator 
tips, and frequency corresponding to the blade passage fre- 
quency. 
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FIG. C.8. CASE Bl; INEVITABLE RADIATION. 
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C.2 Acoustics of a Moving Medium 


The only task that needs to be performed in this section is 
to investigate how the radiation span (-k^ , +kg^ ) along the 
wavenumber k. gets modified due to the o o fact that the 
acoustic medium is moving uniformly with subsonic velocity 
u = (u, , u,, u,), where u,, u,, u, are the velocity components 
in the X, y and z dlreotlins. 


Once again, consider the line monopole source of Pig. C.l, 
with Eqs. (C.l) through (C.5) still applicable. In addition to 
the radial coordinate r of Eq. (C.3), the corresponding radial 
vector r is defined as 


r = (y,z) 


(C.37) 


The far field pressure P(x,r,t) now must satisfy the follow- 
ing field equation, (Morse and Ingard, 196^4), 




1 + 


1 /„ 9 d. 9 d. 9 \ 

0) \ 1 9x ^2 3y ^3 3z / 


p = 0 


(C.38) 


V 


fl 


For Uj = u^ = Uj = 0, Eq. (C.38), reduces to the usual Helmholtz 
equation applicable for a stationary acoustic medium. In analogy 
with Eq. (C.ll), a criterion is needed that the radial wavenumber 
kp must satisfy for no radiation to occur to the far field (l.e.. 
Case 2). However, in analogy with generalization of Eq. (C.3) to 
Eq. (C.37), a radial wavenumber vector kj. is defined as 


= (k. 




(C.39) 


where k^ and kj are the wavenumber components of the radially 
outwards propagating wave. 

n 

Now the important phase aspect of the far field pressure 
P(x,r,t) (for a given value of r in the farfleld) is given by 
the correspondingly generalized version of Eq. (C.6) 

P(x,r,t) = constant x e ° ° , (C.4Q) 

If 


C-l6 


n 



where 


k *r=ky + kz 
r — 2 3 


(C.4D 


The following definitions are required: 


so that 


where 


k = 1 k_ I = (k^ + k^ ) 


u_ = |u„| = (u? + uM 


^2 = k^slnuj^ , 


k = k cosa, , 
3 r k * 


u = u sina , 
2 r u * 


U3 = u^cosa^ , 


k *u =ku +ku = k u^cosa , 
r r 22 33 r r 


(C.42) 

(C.43) 

(C.44) 

(C.45) 

(C.46) 

(C.47) 

(C.48) 


a = (ttj^ - a^) . (C. 49 ) 

Thus, a is the angle between the radial wavenumber vector k^ (or 
the radial location vector (y,z) of observation in the far— 
field) and the radial flow vector Uj,. 

Now, substituting Eq. (C.40) into Eq. (C. 38 ) gives the follow- 
ing relation, 


C -17 



N 



■'1 


[1 - 


0) 


(u k 
1 o 


+ u^k^cosa)]^ 


0 


(C.50) 


I 


which can be rewritten as a quadratic equation in k as follows 

r 


Ak^ + Bk + C = 0 
r r 


(C.51) 


where 


A = 


k 2 
^o 

u^ cos^a - 1 = cos^a - 1 

0) 2 r r 


(C,52) 


B = 


k^2 2 

a a 

2 : u k^u^ cosa - 2 — - 

,, 2 1 o r 0 ) 

w. o 


u^cosa 

r 


2(m m^k^cosa -km cosa) 
1 r o r 

0 


(C.53) 


and 


C = 


(-k^ + k ^ 
° ^o 


- 2 


k ' 
a 


k 


a 

u k^ + —2 — k 2) 
‘‘Jq 1 o 0) 2 o 


= (-k^ + k^2 - 2k„ m k^ + m^k 

o a^ a^ j o , o ' 


(C.54) 


[T 


f 


In the above equations appropriate Mach numbers are introduced 

by division of velocities by the sound speed c = w /k 

o a„ 




[ 
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For the stationary acoustic medium, the condition on radial 
wavenumber magnitude kj., for no radiation to occur (Case 2), was 
that kj, be imaginary (‘see Eq. C.ll). In analogy with that re- 
quirement for no radiation to occur, it is required that kj, of 
Eq. (C. 5 I) be complex (with positive imaginary part). That will 
happen if and only if 


\ 


AC - B^/i| > 0 


(C.55) 


Now, the left hand side of Eq. (C.55) does not contain k^, 
but is 'a quadratic form in kQ, the wavenumber of the source 
wave. Thus, Eq. (C.55) can be rewritten as 


where 


D 

E 

P 


DK^ + Ek^ + P > 0 , 

0 0 ^ 

(C. 56 ) 

1 - m? - m^ cos^a 
1 r 

(C.57) 


(C.58) 

“'1 

(C.59) 


b 


The minimum value of D occurs for a = 0 or tt (i.e., when kj, 
and Up are coincident or oppositely directed. This minimum — 
— value Is given by 


D 


mlnm 




1 - m^ > 0 , 


where m is total flow Mach number. Since and therefore D 

is always positive, the left hand side of Eq. (C. 56 ) is 

positive for large jk^l (i.e., for k^ -*• ±“), being dominated by 
the first term Dki. This behavior, incidently, is consistent 
with the inequality expressed by Eq. (C.56). Recalling that 
the inequality of Eq. (C.56) is the condition on ko for no 
radiation to occur (i.e.. Case 2), the radiation will, in fact, 
take place for a range of wavenumbers k^ of relatively small 
magnitude. This range, the radiation span along axial wave- 
number kj , is determined by the two real roots k„ and k^ of 
the quadratic form of Eq. (C. 56 ) 0 + o- 


C -19 


] 



(C.60) 


If 



-m ± ( 1-m^cos ^a) 
k — i ^ 

( 1-m^-m^cos ^a) 

1 r 


In analogy with Pig. C.2, this radiation span is shown as a 
shaded strip in Pig. C.9. The center O' of the span is shifted 
to the left by the amount k^ mj/D (with D given by Eq. (C.57) 
from the origin 0 (k^=0). o This shift resulting from first 
(common) term on the right hand side of Eq. (C.60), is interpreted 
as a Galilean shift. The equal intervals (k , O') and O', k^ ), 
resulting from the second terms on the right “o- hand side of O’*" 
Eq. (C.60) are Interpreted as Lorentz half-spans. 

Also shown in Pig. C.9 is q(ki,t) for a phase wave whose phase 
velocity cq is supersonic in the fixed coordinate system (or 
equivalently whose kQ is less than kg^ ), yet since kQ lies out- 
side the radiation span, the ® particular phase wave 

Illustrated will nOt radiate to the farfleld. 

Reformulation of the condition of no radiation in terms of 
Mach numbers can be done along exactly similar lines as done in 
Eqs. (C.19), (C.20) and (C.21). Thus, the upper and lower permis- 
sible Mach numbers m^ and m^^ are given by 

( 1-m^-m^cos ^a) 

„ _ 1 r 

-m -(1-m^cos^a) 

1 r 

(1-m^-m^cos^a) 

_ _ 1 r 

m = 

^ -m^+ ( 1-m^cos ^a) 


and for no radiation to occur, the following must be satisfied: 


T7 = m -(1-m^cos^a) 

/2 1 r 


(C.6l) 




= m +(l-m^cos^a) 


Va 


(C.62) 


m. 


< m < m 


u 


(C.63) 


cr 
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FIG. C.9. 


SKETCH OF RADIATION SPAN ALONG WAVENUMBER 
FOR A MOVING ACOUSTIC MEDIUM. 



Extension similar to that from Figs. C.2 to C.3 can also be 
easily performed for the present case; as a result of the medium 
motion, the acoustic "cones” of radiation in the w-k plane will 
be asymmetrical about the cu axis. 

Finally, the entire discussion of discrete arrays in Sec. C.2 
can be applied here with the newly defined radiation span. 


C.3 An Estimate Of Overall Power Radiated From The Stator 

Figure C.IO shows the wake velocity deficits as viewed in 
time at one SBLE near the tip. There are f such deficits per 
second, where f = 8600 Hz is the rotor blaSe passage frequency. 
Each individual wake deficit, v(t) has an approximately Gaussian 
shape around its peak deficit value v^, thus 

v(t) :: v^ g-tV2T^ ^ (0.64) 


where the "standard deviation" T, and the maximum deficit v are 
estimated to be: ° 


T = 9.6 X 10~* sec. 


(C.65) 


V =44 m/sec (l44 ft/sec) 
o 


( 0 . 66 ) 


The maximum deficit, v corresponds to 10° change in angle of 
attack. The time interval t between consecutive deficits is 
given by 


T = ^ = l.l6 X 10 ** sec. 
^ r 


(C.67) 


Note that x is about an order of magnitude greater than T, 
in other words the wake deficits are narrow in time when compared 
to their rate of arrival. 


The above data regarding the wake velocity deficits was 
developed from Kemp and Sears (Ref. 18). The description of the 
wake velocity deficit in a spatial coordinate, x, can be obtained 
by assuming that the wakes arrive at the SBLE tips as (locally) 
frozen spatial patterns, being convected along with the local gas 
speed V, where 

V= (m ^ c = 195 m/s (64l (C.68) 

ft/sec) 
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li 



FIG. C.IO. SKETCH OF TIME HISTORY OF WAKE VELOCITY 
DEFICITS AS THEY IMPINGE ON A SINGLE 
SBLE TIP. 
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and being given by Eqs. 2,1 and 22 . Thus, the spatial 

picture of a wake' velocity deficit is obtained by the trans- 
formation. 


X = Ut . 


(C.69) 


For estimating the overall acoustic power radiated from 
the stator it is convenient, as shown below, to integrate the 
results in the time domain. However, in order to get a qualita- 
tive understanding of the situation in the frequency domain we 
discuss briefly the Fourier transform v(u)) of v(t) of Eq. (C.6^). 


v(u) 


' ^ 1 '' 


(t)e^“‘ dt 


(C.70) 


Tv 


o ^-a)^T2/2 


(27t) 


Va 


(C.71) 


The Fourier transform v’ (o)) of the sequence of pulses of 
Fig. C.IO may then be written as 


V ' (o)) 


1 

2it 


v’(t)e^“^ dt 


1 

2it 


+ 00 

£ v(t-T.1) e^*^^ dt 

J=_oo 


(C.72) 


+ 00 

= v(w) 0 )^ 2 6(o)-mo^) 

n=-oo 


(C.73) 
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where u is the blade passage frequency in radians/sec (see 
Eqs. 17 and lo 


0). = 2Trf 
o r 



(C.74) 


Thus, as expected, the frequency content of the rotor wake 
velocity deficits v'(t) consists of the various rotor harmonics 
n. Since y(o)) does not decay appreciably with increasing 
frequency (the "standard deviation" of v(w) is 1/T) the higher 
rotor harmonics at n = ±2, ±3, etc. have nearly the same strength 
or amplitude as the fundamental harmonic at n = ±1. 


Now, reverting back to the time domain analysis, the 
fluctuating lift A(t) generated at the leading edge of a SBLE 
tip from impingement of one wake deficit v(t) is given by 


^'(t) = I v(t’) h(t-t') dt ' , (C.75) 

0 

where h(t) is the unit Impulse response function derived from 
Kvissner's function [Blsplinghoff et al. ^ Ref. 19]. Since the 
essential uncancelled fluctuating lift is restricted to a 
relatively small segment of the SBLE span near the tip, use 
of Kvissner's function, valid for low aspect ratio, is readily 
justified for the present calculation. Since Kvissner's function 
gives the lift due to a sharp-edged gust (i.e., due to a gust 
that is a unit step function), the unit impulse response function 
h(t) is obtained by differentiating the Kvissner's function. 
h(t), so obtained, is given by 


h(t) = 


0.13 



1 


e 


5 


(C.76) 


where t. is the time taken by the gust to travel the (swept) 
semichord b. 


= b/U 5 1.33 X lO”** sec. (C.77) 


The above estimate of tk is based on b ~0.03^ m (0.11 ft), and U of 
Eq. C.68. The lift coefficient C, is given by 


Cl ~ X » (C. 78 ) 


where p is the medium density ( 2.4 x 10~* lb-sec ^/ft“ 

= 1.24 Kg/m^) and Xj, is the acoustic wavelength at blade passage fre- 
quency fj, (Eq. 26 ). Note that in Eq. C. 78 , Xp/4 denotes a rough 
estimate of the SBLE span near the tip from which the un- 
cancelled fluctuating lift is estimated to radiate. Now, this 
choice of X^/4 is suitable (only) for the rotor fundamental 
harmonic at frequency u . For the higher rotor harmonics of 
Eq. C. 73, correspondingly smaller spanwise length scales would 
be more appropriate. However, in the time domain analysis that 
is being pursued, the choice of X /4 in Eq. C .78 is taken to 
apply to all the rotor harmonics, therefore the resulting 
estimate of the overall (l.e., frequency-integrated) radiated 
power is liable to be conservative. 

Substituting Eqs. (C.64) and (C.76) into Eq. (C.75), enables 
calculation of fluctuating lift il(t) at a single SBLE tip due to 
the impingement of a single wake velocity deficit v(t). Since 
the minimum time constant of h(t) is much larger than the time 
constant or "standard” deviation” T of v(t) [compare Eqs. (C.65) 
and (C.77)], for evaluating il(t) from Eq. (C.75), one can justifi- 
ably approximate v(t) of Eq. (C.64) as follows, 

v(t) = V (2 tt )^^2 T6(t) . (C.79) 

0 
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The constant (2tj) T In Eq. (C.79) Is Introduced so as to make 
the total "area" (in other words, the integral /v(t) dt the 
same for Eqs. (C.64) and (C.79). From Eq. (C.70) and (C.7D, note 
that this area is equal to 2ttv(£o)i Substituting Eq. (C.79) 

into Eq. (C.75), gives I 


t(t) = Vj(27t) T h(t) . (C.80) 

I 

Prom the point of view of generation of steady lift, the 
stator blade chord is expected to be oriented parallel to the 
flow velocity U at its leading edge, so that a zero mean angle 
of attack is ensured. Hence, the wake-deficit-induced fluctuating 
lift of Eq. C .80 is oriented normal to the flow velocity U. The 
acoustic intensity I(t) radiated by this "transverse" dipole 
(i.e., the direction of fluctuating lift is normal to flow) is 
given by [Lighthill ,Ref . Ig; Morse and Ingard, Ref. l7], 


Kt) = [2A(t)]2 L_ G (jn) 

12p7TC^ ^ 


(C. 81) 


where £,(t) - d/dt £(t) and G 2 (m) is a function of the flow Mach 
number m = U/c , 




2 

m^d-m^) 



1+m 

1-m 


(C.82) 


The factor 2 appearing with i{t) in Eq. C .81 accounts for the 
baffling effect due to the presence of the duct wall (assumed 
to be acoustically rigid) enclosing the SBLE tip. 


The radiated acoustic energy E, associated with the intensity 
I(t) of Eq. C .81 is given by 


E 


tCO 

iCt) dt. 

0 


CC.83) 


M 


V 
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N 


] 


From Eqs . C .76, C.80 and C.8l, we note that the only time depen- 
dent factor of I(t) Involves h(t) of Eq. C. 76 , hence, the Integral 
to be evaluated is 


h(t)^dt = C 2 , 


\ 


thus, substituting Eqs. C. 80 , C. 8 I and C. 84 into Eq. C .83 vie get 


E 





pc 


(v TCt 
0 


0.133 



(C. 85 ) 


Now, E is the acoustic energy radiated from a single SBLE 
tip due to impingement of a single wake velocity deficit v(t). 
Hence, the acoustic power radiated from a single SBLE tip is Ef , 
where f is the rate of impingement of wake deficits on the ^ 
SBLE tip. 


1 Next, assume that the V individual SBLE tips radiate 

< more or less incoherently (an assumption particularly valid 

for higher rotor harmonics n of Eq. C. 73). Hence, the power 

radiated from the V stator tip sources is Ef V. 

r 

Finally, even though the calculation for the power radiated 
from the stator leading edge sources at the hub is not carried 
out separately, because of closer circumferential spacing between 
these hub sources, the total power radiated from the stator hub 
) is likely to be considerably less than that from the stator tip. 

The total power II radiated from the stator is conservatively 
estimated to be given by 


) 


n = 2Ef V 
r 


(Cv86) 
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Substituting In Eqs . (C.85) and (C. 86 ) the numerical 
quoted for various parameters (with V = 59 ) gives 


n - 130.5 dB re 10 watt . 


values 

(C. 87 ) 
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APPENDIX D: NOTES ON EMPIRICAL CALCULATION OF FAN NOISE LEVELS 


As mentioned In Section 7> all components of the rotor and 
stator noise spectrum could not be calculated from basic con- 
siderations. The previous Appendix gives a noise level cal- 
culation for the residual stator noise sources. In the Interest 
of determining what reduction In levels the swept rotor might 
be expected to cause, Burdsall's empirical correlation (Ref. 

19) » was exercised for both the actual fan model and a "full 
scale" counterpart. The parameters required In Burdsall’s 
routine are given In Table D-1. Figure D-1 summarizes the 
narrow band power levels and spectra for the various components, 
(SPL arbitrarily computed at I 50 ft., 60° from rotor axis) 
showing the predominance of MPT'S. Of course, the details of 
the MPT spectrum vary from fan to fan due to their origin In 
manufacturing tolerances. Pig. D-2 shows a typical comparison 
of Burdsall’s prediction with measured data. Indicating a fairly 
large fluctuation In harmonic levels around the mean line of 
the prediction. In Fig. D-1, It Is clear that according to this 
scheme, elimination of MPT’s would reduce the tone levels 
considerably. However, note that In Pig. D-1, the line Is an 
envelope of discrete frequency levels at various multiples of 
rotation rate while the broadband spectrum is contlnous. Thus, 
integration Into constant percentage bandwldths, and into overall 
levels will lead to the MPT and broadband levels beipg very nearly 
identical . 

As a final point, it is Interesting to note that the power 
level of the BPP tone (non-MPT noise) is ~10 dB above the 
predicted level for the swept stator as described in Appendix 

C. 
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APPENDIX E. ALGORITHM FOR DERIVATION OF STATOR LEADING EDGE TRACE i 

VELOCITY IN STATOR FIXED COORDINATES 


The geometry of a rotor wake as it reaches the stator is 
given in Pig. E-1. Refer to Pig. 9b in the text for 3-dimensional 
representation. The following four steps give the rotor wake shape 
and local trace velocity for both an unswept stator (or at the 
inlet plane of a swept stator), and for swept stators. Aerodynamic 
reaction on the rotor path by the stator is not taken into account. 



FIG. E-1. GEOMETRY FOR CALCULATION OF ROTOR WAKE SHAPE AND TRACE 
SPEED ON STATOR VANES. 


1) Unswept, Unskewed Stator 

The skew of the rotor wake at stator plane is a(r). The 
local angle between the wake and the radial direction may be 
derived from: 


tan 6 = r ^°‘/6r 

w 


z 


[ 
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The trace velocity In the radial and axial directions is 
respectively 


\ 


(r) 

Br 


(wr)/tan ^(r) 


and 



0 


where : 


V = the trace velocity in blade fixed coor- 

■^Bj^ y. dlnates for radial and axial directions 

* respectively. 


2) Add Sweep Pg 

The rotor wake in the conical surface of the swept leading 
edge is changed, as follows. 


(r) = 


■B 


R 


wr 


tan6 ’ (r ) 
w 


where tan6 * = r 4^' 

w ^ 


and 


V,j,_(r) = tanug(r) where UrCt) = local blade sweep angle 


^B 




B' 
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a'(r) = a(r) + Ax 


> 


where Ax Is the downstream displacement of the leading edge caused 
by sweep. 


tane ’ 

w 


(r) = r- 


g(r) 


5g (r) 
+ Ax (Sx 
6r 


= r <Sa(r) 
6r 


+ r Ax 


<S^ot(r) 6Ax 6g(r) 

6x6r 6r 6x 


= tan0 + 
w 


6tanew . _ 6g(r) 

6 X ^ ~53T" 


tanug 


tan0'^(r) = tan0^ + r tanu^ 


3) Trace Velocity for Swept Stator in Stator-Fixed Coordinates 


The radial component of trace velocity Is 


(r) 

Br 


wr 


tan0 


* + 


w 


ZT 


tany 


B 


* <5 tan 

where tane* (r,x) = tan0^ + Ax — 


The axial component Is 


9 


m (r) 


wr tan0g 
tanej . 


9 


tanUg 



where 0* is determined from cross-plots of the wake path in the 

r,a plane at various axial locations, and 6a/6r = the local wake 
helix pitch angle determined from wake crossplots in the meri- 
dional plane. 


4) Transformation of Trace Velocity Amplitude From Stator- 
Fixed Coordinates to Gas-Fixed Coordinates 



where Vq Is the gas velocity in the axial and circumferential 

x,c 

directions, respectively. 
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LIST OF SYMBOLS 
axial 

inlet area (to rotor passage) 

geometric throat area 

area to choke the flow 

normalized distance 

blade semichord 

number of blades 

chord length; sound speed 

phase or trace velocity; sound speed 

center of gravity 

lift coefficient 

diffusion factor 

double circular area (blade profile) 

distance from blade section c.g. to pressure 
surface; circumferential spacing between 
adjacent stator tips 

swept leading edge element 

acoustic energy (radiated from a single SBLE tip) 
Mach factor ( = 1/M ) ; frequency 

blade passage frequency of rotor blade 


S-1 









LIST OF SYMBOLS (Cont.) 

■ 

AFj . 

= 

centrifugal force at center of a blade volume 
element located at j 


Gj (m) 

=z 

function of flow Mach number 


h 

= 

unit Impulse response function 

II 

I 

= 

acoustic Intensity 


j 

= 

source number 


k 

= 

wavenumber; constant 


k 

r 

= 

radial wavenumber 


L 

= 

harmonic 


1 

= 

distance from section c.g. to leading edge; 
fluctuating lift 


LCP 

= 

low cycle fatigue 


LE 

= 

leading edge 

1 

mn 

= 

rotor/stator Interaction harmonic 


m 

— 

circumferential mode number; component of 
Mach number; function In deviation angle 
formula 


M 

.= 

Mach number; moment 


"w 

= 

relative flow Mach number component 

m 

"w. 

= 

relative Inlet Mach number 


M » 

Wi 

= 

Mach number required for a subsonic edge to 
achieve sonic ( = M ) 

Wi 

s: 

M 


component of Mach number which Is always normal 
to the leading edge ( = 1 for sonic LE; < 1 
for subsonic LE 
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LIST OF SYMBOLS (Cont.) 


M. 


w 


'L 


Mach number required to make a subsonic LE 
a sonic LE 


am 


IJ 


N 


n 


P 

P or p 
P 

P/P 

P 

PNL 

q 

A 

R 

«c 

r 


s 

SAP 

SCP 

SBLE 

SR 


moment of j force about c.g. of section i 
rotation rate 

shape parameter for polynomial blade forms 
harmonic number 

total pressure; static pressur 
far field acoustic pressure 
location of leading edge point 
pressure ratio 

static pressure after normal shock 
Perceived Noise Level 
source strength 

monopole source strength per unit length 

distance from origin 

radius of curvature of streamline 

radial distance 

circumferential blade spacing 

Structural Analysis Program 

Stress Concentration Factor 

Stator Blade Leading Edge 

sweep reversal 
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LIST OF SYMBOLS (Cont.) 


t 

t (x) 
T 

TE 


U or u 

V 

Vo 



W 

Was 

w/w 

w 

Ax 


X 

z, Zl 


thickness; time 
thickness distribution 
standard deviation (time) 
trailing edge 
air velocity; wheel speed 
velocity 

peak velocity defect 

number of stator vanes; velocity 

volume element of a blade 

mass flow rate; velocity; inlet relative 
velocity 

average velocity 
specific mass flow 

flow deceleration rate 
relative velocity 

downstream displacement of SBLE caused by 
sweep 

linear distance 

axial coordinate of leading edge 



LIST OF SYMBOLS (Cont.) 


GREEK 


angle of attack; angle between radial wave 

vector k and radial flow vector u ; local 
r r 

Mach angle; wake displacement angle from 

radial 

stator Inlet angle 

relative flow angle; exit angles of flow 
setting angle; ratio of specific heats 
flow deviation angles; Dirac delta function 
stator deviation angle 

slope of the relative flow velocity; slope 
between lines connecting section LE and CG 
and a line connecting LE with lower surfaces 
coordinate at mid-chord (e^ = sin“^Vj,/wi) 

centerline-projected displacement of the 
c.g.'s of an airfoil section at r. relative 
to one at r^^ ^ 

acoustic wavelength 

acoustic wavelength at blade passage frequency 
polytropic state efficiency 
circumferential angle 

angle between radius and rotor wake centerline 
unswept stator 

angle between radius and rotor wake centerline 
swept stator 









LIST OF SYMBOLS (Cont.) 

V 

®L 

= 

angular abscissa of leading edge point 


y 

= 

Mach cone angle; radial order of acoustic 
modes 



= 

stator sweep angle 

IT 

y” 

= 

projection of the Mach cone angle on the 
w-r plane 

1 & 

V 

s 

hub-to-tip ratio; lateral sweep angle; 
section thickness ratio (t /c); relative 

thickness 


n 

= 

acoustic power 


p 

= 

density 


a 

= 

sweep angle; stress; cascade solidity 


T 

= 

shear stress; LE and TE thickness factors; 
time Interval between successive events 


"b 

= 

time for a gust to travel the swept semi- 
chord (b) 

r 

<P 

— 

local camber angle; angle between leading 
and trailing edge along the unwrapped conical 
surface 


$ 

= 

mean camber angle 


0) 

= 

radian frequency; wheel rotation speed 

a 

0) 

r 

= 

radian frequency of blade passages 




shaft rotation frequency 

u 
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LIST OF SYMBOLS (Cont.) 




SUBSCRIPTS 

A 


axial; along blade leading edge 

a 

= 

acoustic 

ax 

s 

axial 

B 

= 

blade-fixed coordinates 

c 

= 

circumferential 

crit 

= 

critical 

eg 

= 

center of gravity 

D 

= 

defect 

G 

= 

gas-fixed coordinates 

g 

= 

geometric 

i 

= 

component parallel to blade array 

i, J 

= 

indices of spatial coordinates 

L 

= 

normal to leading edge; leading edge 

1 

=: 

lower 

M 

= 

mainstream 

m 

= 

meridional component 

mi run 

= 

minimum 

0 

= 

trace speed 

R 

s 

radial 


LIST OF SYMBOLS (Cont.) 


radial; component normal to blade array 

tangential; wake trace 

tip 

upper; tangential 
tangential 

chordwise distance from LE; direction normal 
to z-axls 

Cartesian coordinate normal to z-axis 
axial 

freestream relative 
X, y, and z directions 
rotor inlet station 
rotor exit station 
stator inlet station 


stator exit station 


